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Final report on grant NAG5-3039 


Grant NAG5-3039 includes the proposals: 

1. The Complex X-ray Spectrum of Quasars with Soft Excess 

2. Ionized Absorbers in AGN 

3. Very Steap Soft X-ray Quasars 

The activity concerning the three proposal ended. 

A paper presenting the results of the ASCA observation of NGC985 has been published 
on The Astrophysical Journal (Nicastro, Fiore k Matt 1999, ApJ, 517, 108) and on the 
proceedings of the Conference “The Active X-ray sky” (Nicastro, Fiore, Brandt k Reynolds, 
1998, Nuclear Physics B, 69/1-3 (1998) 501. This observation showed the importance of 
resonant absorption in modifying the 0.5-1. 5 keV spectrum of this source. Detailed physical 
models have been produced. 

The ASCA observations of the two Soft X-ray Quasars PG1244+026 and NAB0205+024 
have been published on the Montly Notices of the Royal Astronomy Society (Fiore et al. 1998, 
MNRAS, 298, 103). A spectral feature at about 1 keV has been discovered in PG1244+026, 
similar to those in other steep X-ray spectrum Seyfert 1 galaxies and quasars. The line 
could result from reflection from an highly ionized disk, in line with the view that steep 
X-ray spectrum AGN are emitting close to their Eddington luminosity. 

The papers quoted above are attached to this report. 


Best Regards, 
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ABSTRACT 

Steep soft X-ray (0.1-2 keV) quasars share several unusual properties: narrow Balmer lines, 
strong Fen emission, large and fast X-ray variability, and a rather steep 2-10 keV spectrum. 
These intriguing objects have been suggested to be the analogues of Galactic black hole 
candidates in the high, soft state. We present here results from ASCA observations for two of 
these quasars: NAB 0205 4- 024 and PG 1244 + 026. 

Both objects show similar variations (factor of ~2 in 10 ks), despite a factor of 
—10 difference in the 0.5-10 keV luminosity (7.3 x 10 43 erg s -1 for PG 1244 + 026 and 
6.4 xlO 44 erg $~ l for NAB 0205 4-024, assuming isotropic emission, H 0 = 50.0 and 
q 0 = 0.0). 

The X-ray continuum of the two quasars flattens by 0.5-1 going from the 0.1-2 keV band 
towards higher energies, strengthening recent results on another half-dozen steep soft X-ray 
active galactic nuclei. 

PG 1244 4- 026 shows a significant feature in the ‘1-keV’ region, which can be described 
either as a broad emission line centred at 0.95 keV (quasar frame) or as edge or line absorption 
at 1.17 (1.22) keV. The line emission could be a result of reflection from a highly ionized 
accretion disc, in line with the view that steep soft X-ray quasars are emitting close to the 
Eddington luminosity. Photoelectric edge absorption or resonant line absorption could be 
produced by gas outflowing at a large velocity (0.3-0.6 c). 

Key words: line: formation - line: identification - galaxies: individual: NAB 0205 4- 024 - 
galaxies: individual: PG 1244 -h 026 - galaxies: Seyfert - X-rays: galaxies. 

(ii) strong Fe u emission (Laor et al. 1994, 1997; Lawrence et al. 
1997); 

(iii) rapid, large-amplitude variability - factor of 2-50 on time- 
scales from minutes to months (Boiler et al. 1 996; Brandt, Pounds & 
Fink 1995; Otani, Kii & Miya 1996; Boiler et al. 1997); 

(iv) somewhat steep hard X-ray spectra, 2 > a 2 -iok e v > 0.6 
(Pounds, Done & Osborne 1995; Brandt, Mathur & Elvis 1997). 

Pounds et al. (1995) suggest the latter to be a close physical 
analogy with the X-ray power law produced by Comptonization in a 
hot accretion disc corona in Galactic black hole candidates (BHCs) 
in their ‘soft, high’ state. This is not the only analogy between BHCs 
and steep X-ray spectrum quasars. Laor et al. (1994, 1997) 
explained the correlation with H/3 full width at half-maximum 
(FWHM) to be the result of the larger size of a virialized broad 
emission-line region for an active galactic nucleus (AGN) in a high 
L/Lzm state. Ebisawa (1991) found that, while the soft component 

© 1998 RAS 


1 INTRODUCTION 

The ROSAT Position Sensitive Proportional Counter (PSPC) has 
found a large spread in the energy spectral indices of low-z 
quasars 1 : 0.5 < ao j_ 2 keV) < 3.5. In about 10 per cent of cases 
<*ai- 2 kev ~ 2 (e.g. Laor et al. 1994, 1997; Walter & Fink 1993; 
Fiore et al. 1994). The large spread in a 01 _ 2ktV favoured the 
discovery of its correlation with other properties. In fact, the steep 
soft X-ray quasars have then been realized to share a cluster of 
unusual properties: 

(i) narrow Balmer lines 2 (Laor et al. 1994, 1997; Boiler, Brandt 
& Fink 1996); 

1 We use ‘quasars’ to describe broad line emission objects, regardless of 
luminosity. 

2 The permitted lines have FWHM S 2000 tons" 1 , yet still are clearly 
broader than the forbidden lines. 
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of six BHC observed by Ginga is roughly stable on time-scales of I 
day or less, the hard component exhibits large variations down to 
millisecond time-scales. These time-scales translate to 10 4 yr and 
0. 1 d for quasars, if they scale with the mass of the compact object 
The soft component of a BHC extends up to — 10 keV in a BHC in 
‘soft, high’ states, and it is often associated with optically thick 
emission from an accretion disc. If this is the case, the temperature 
should scale with the mass of the compact object as A/bh M * and the 
above energy translates to 0. 1 -0.4 keV for quasars. The rapid large- 
amplitude variability shown by a few narrow-line Seyfert 1 galaxies 
(NLSy Is) at about 1 keVon time-scales of hours to days (Otani et al. 
1996; Brandt et al. 1995; Boiler et al. 1997) can then be analogous 
to the above BHC hard component flickering. 

A steep X-ray spectrum quasar with 10-100 times the luminosity 
of NLSy 1 s should be larger and so should vary no more rapidly than 
several days. Instead Fiore et al. (1998) find that steep X-ray 
spectrum PG quasars commonly vary by a factor of 2 in one day. 
Variability therefore seems to be correlated with X-ray spectral 
slope and Balmer linewidth (and therefore possibly with the 
accretion rate) rather than with the luminosity. 

Evidence for spectral features in the * 1 -keV’ region in many steep 
soft X-ray quasars is building up (Turner et al. 1991; Brandt et al. 
1994; Otani et al. 1996; Comastri, Molendi & Ulrich 1995; Leighly 
et al. 1997, in preparation). Instead, ‘normal’ Seyfert 1 galaxies 
(having broad Balmer lines and flatter soft X-ray spectra) usually 
have their strongest absorption features at lower energies (in the 
0.6— 0.9 keV ‘oxygen’ band). An intriguing possibility is that the 
appearance of these features al different energies also depends on 

Ul** i- 

Detailed high-energy X-ray spectra of luminous quasars with 
steep soft X-ray spectra are essential to understand the ‘ narrow - 
broad line’ phenomenon in AGN, in particular whether the peculiar 
X-ray properties depend on optical luminosity, on optical-to-X-ray 
ratio (a ox ), or on their Eddington ratio. To this end we selected 
two bright quasars with oo i_ 2 k e v >2.0 (Fiore et al. 1994) at 
the extreme values of optical luminosity, both with low Galactic 
V H (Table 1) of 1.9 XlO 20 cm' 2 for PG 1244 + 026 and of 
3.0 x 10 20 cm”' for NAB 0205 + 024 (Elvis, Lockman & Wilkes 
1989) and observed them with ASCA (Tanaka, Inoue & Holt 1994). 
We report the results in this paper 

2 OBSERVATION AND DATA REDUCTION 

Table 1 gives the redshift, M Vf the 0.2-2 keV luminosity, ao X , the 
average PSPC count rate and spectral index and the Galactic N H for 


the two quasars. Table 2 gives the ASCA observation log, the SIS 
and GIS exposure times and count rates. 

Both observations were performed in two charge-coupled device 
(CCD) mode with the source at the ‘ICCD mode’ position. Data 
reduction was performed using ftools 3.6. We used ‘bright’ mode 
SIS data, combining low, medium and high bit rates. Conservative 
cleaning criteria were applied (minimum Earth occultation = 7°, 
minimum magnetic rigidity = 6 GeV/c, minimum bright Earth 
angle = 20°, and excluding data collected in the first 32 s after 
the satellite passage in the South Atlantic Anomaly (SAA) and 
through the day -night terminator). Counts, light curves and spectra 
from the two quasars were accumulated in circular regions of 3 and 
4 arcmin radius for SIS and GIS respectively. 

We are interested in the high-energy spectrum of these sources 
and since they are rather faint, and possibly very steep, background 
subtraction plays a crucial role. Background counts were accumu- 
lated from regions surrounding the sources and compared with 
counts accumulated from the same regions from ‘blank-sky’ 
observations. The ‘local’ and *blank-sky’ background counts 
were always within 10 per cent of each other for the four ASCA 
instruments. To obtain the best possible signal-to-noise ratio in the 
background-subtracted spectra, we therefore used the ‘blank-sky’ 
background in our spectral analysis. We extracted background 
spectra from ‘blank-sky’ event files using the same regions as for 
source extraction. The count rates of the two sources become that of 


03 PG 1244+026 
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Figure I. The SIS 0.5-7 keV light curve of PG 1244 + 026. 


Table 1* Steep soft X-ray quasars. 


name z 


*■0.2-2 

a ox 

PSPC a o.l-2 keV ^H(Gan 



10 45 erg s" 1 


cts/s 

10 20 cnrT 2 

PG 1244+026 0.048 

-21.1 

0.14 

1.4 

1.0 2.3±0.1 

1.9* 

NAB0205+024 0.155 

-25.0 

1.8 

1.6 

0.7 2.3±0.1 

3.0* 

*H 0 = 50, q 0 = 0. b Elvis et al. (1989). 






Table 2. ASCA observations. 






name Dates 


Exposure SIS-GIS 

2 SIS count rate 

2 GIS count rate 



ks 


cts/s 

cts/s 

PG 1 2*4+026 1 -3 J u 1 1 996 

37-39 


0442+/-0.004 

0.241 +/-0.004 

NAB0205+02* 18-20 Jan 1996 

50-54 


0.150+/-0.001 

0.153+/-0.002 
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the background at about 7 ke V (observer frame). After background 
subtraction PG 1244 + 026 is observed in the GIS up to 10 keV and 
NAB 0205 + 024 up to 8 keV (9.3-keV quasar frame), both at the 
> 3a level. 

Spectral fits were made separately to the spectra from the four 
ASCA instruments and to the spectra obtained combining together 
the data from the two SIS and GIS detectors. The results were 
consistent with each other. In the following we present the results 
obtained following the second approach. In some cases \ are 
smaller than 1. This is because of the prescription adopted in 
adding the spectra, for the propagation of the errors (the Gehrels 
1986 algorithm: error — 1.0+ SQRT(N + 0.75)). Spectra were 
always rebinned following two criteria: (a) to sample the energy 
resolution of the detectors with four channels at all energies 


NAB0205+024 
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Figure Z The SIS 0.5-7 keV light curve of NAB 0205 + 024, 


ASCA observations of two X-ray quasars 105 

where possible, and (b) to obtain at least 20 counts per energy 
channel. In the spectral fits we limited ourselves to the 0.6-10 keV 
energy band, to minimize the systematic effect caused by the 
uncertainty in the SIS calibration below 0.6 keV. Af H was always 
constrained to be greater than or equal to the Galactic value along 
the line of sight 

In all cases the quoted errors represent the 90 per cent confidence 
interval for one interesting parameter. 

3 VARIABILITY 

ASCA observed the two quasars for a total elapsed time of 105 ks 
(PG 1244 + 026) and 1 19 ks (NAB 0205 + 024). We can therefore 
study the variability of these sources on time-scales from a few 
hundred seconds to about 1.5 d. 

Figs 1 and 2 show the SIS light curves of the two sources. 
Following Ptak et ah (1994) the values were computed using 
variable bin sizes corresponding to good time intervals longer 
than 200 s. 

Significant variations are evident in both light curves. A factor 
of 2 variability, both down and up, in about 10 ks is present at the 
beginning of the light curve of PG 1 244 + 026. A drop of a factor of 
—2 is present at the beginning of the light curve of NAB 0205 + 024 
on a ~ 15-ks time-scale. In this paper we limit ourselves to pointing 
out that roughly similar variability is observed in two sources that 
differ in luminosity by a factor of ~~10. Significant, although rather 
small, spectra! variability is also present in the ASCA observations 
of the two sources. A systematic analysis of the variability in 
different energy bands and of the spectral variability is in progress 
and will be presented, together with a similar analysis on a sample 
of about 20 Seyfert 1 galaxies and quasars observed by ASCA and 
BeppoSAX , in a paper in preparation. In the following sections we 
present the average properties of the spectra. We anticipate that the 
spectral variability will not modify the results presented here. 



Figure 3. The SIS + GIS spectra of PG 1 244 + 026 fitted with a simple absorbed power-law model. The lower panel shows the ratio between the data and the 
best-fitting model. 
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4 SPECTRAL ANALYSIS 

Figs 3 and 4 show the SIS 4* GIS spectra of PG 1244 + 026 and 
NAB 0205 + 024 fitted with a simple power law absorbed at low 
energy by a column of cold gas equal to or higher than the Galactic 
column along the line of sight. Table 3 gives the best-fitting 
parameters and x". It is dear that this simple model is inadequate 
to describe the 0.6-10 keV spectrum of both quasars A hard tail, 
larger than the 5-10 per cent systematic uncertainties at these 
energies (Gendreau & Yaqoob 1997), is evident in both cases. In the 
spectrum of PG 1244 + 026 there is also a significant excess with 


respect to the model at about 1 keV. The feature is visible in both 
SIS and GIS detectors. We discuss these two findings in turn in the 
next two subsections. 

4.1 The 0.6-10 keV continuum 

The ASCA SIS and GIS responses are peaked at 1-2 keV and 
decrease sharply at higher energies. Therefore the ASCA ‘2-10 
keV' slopes are strongly biased towards the lowest energy bound- 
ary. This means that some caution should be used when comparing 
ASCA ‘2-10 keV' slopes with those of experiments whose 


ASCA SIS - GIS 


NAB0205+024 



Figure 4. The SIS + GIS spectra of NAB 0205 + 024 fitted with a simple absorbed power-law model. The lower panel shows the ratio between the data and the 
best-fitting model. 

Thble 3. PG 1 244 4- 026 and NAB 0205 + 024 spectral fits. 


model 


a £ or T b 

a H or Clf2n 

^.* orAC 

X 2 <dof) 

PG 1244+026 
PL 0.6-10 keV 

1. 9+0.6 

1.67±0.04 



177.9 (148) 

PL 2-10 keV 

1.9FIXED 

1.35±0.12 

- 

- 

48.4 (82) 

PL 3-10 keV 

1.9FIXED 

1.03±0.30 

- 

- 

29.1 (58) 

PL 4-10 keV 

1.9F1XED 

0.67±0.55 

- 

- 

14.4 (38) 

Broken PL 0.6-10 keV 

3.6!?; 7 6 

180 Iq'I 

io6 m 

2.8±0.6 

149.1 (145) 

PL+Raym 0.6-10 keV 

1. 9+0.7 

0.88 %g 6 

1.54+0.06 

>0.5 

135.3 (145) 

PL+BB 0.6-10 keV 

5 . 3 :!.? 

0.16±0.03 

1.4010.15 

- 

139.2 (145) 

PL+Comp.Refl 0.6-10 keV 

2.6^7 3 

2.78±0.06 

>8 

- 

150.5 (146) 

NAB0205+024 
PL 0.6-10 keV 

3.0+0. 1 

1.38±0.02 



126.2 (147) 

PL 2-10 keV 

3.0FIXED 

1.09±0.10 

- 

- 

42.4 (86) 

PL 3-10 keV 

3.QF1XED 

Q.94±0.20 

- 

- 

23.3 (58) 

PL 4-10 keV 

3.0FIXED 

0.62±0.45 

- 

- 

9.5 (38) 

Broken PL 0.6-10 keV 

3. 0+1. 4 

1-50 ig;? 

0.98±0.10 

2.4±0.5 

85.2 (145) 

PL+Raym 0.6-10 keV 

3.0+7. 0 

0.49±0.05 

1.02±0.25 

<0.01 

82.8(144) 

PL+BB 0.6-10 keV 

3. 0+0.2 

0.16+0.02 

1 . 1 5±0.07 

- 

86.2 (145) 

PL+Comp.Refl 0.6-10 keV 

3.0+0.2 

1.44±0.04 

>4.2 

- 

95.2 (146) 


* In 10 20 cm 2 . b In keV. c Metal abundances. 
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responses peak around 6 keV, such as EX O SAT, Ginga, BeppoSAX 
and XTE. To address this, we fitted the SIS 4* GIS spectra of PG 
1244 + 026 and NAB 0205 4 024 with a simple power law in the 
observed 0.6-10, 2-10, 3-10 and 4-10 keV ranges. The results are 
given in Table 3 and shown in Fig. 5. 

In the last three series of fits was fixed to the Galactic value. 
The best-fitting slopes flatten by 0.5-1 going from the low energy to 
the higher energy range. (We note that for NAB 0205 4 024 the 
quasar redshift implies that these slopes refer to slightly harder 
energy ranges: 10 keV in the observer frame corresponds to 1 1.6 
keV in the quasar frame). 



l0 g E Mm 


Figure 5. The PG 1244 4 026 (filled circles) and NAB 0205 4 024 (open 
triangles) best-fitting power-law index in the observed 0.6-10, 2-10, 3-10 
and 4-10 keV ranges as a function of the range minimum energy. The 
NAB 0205 + 024 points have been shifted from real energy for the sake of 
clarity. The PSPC 0.1-2 keV points are also shown for comparison. 
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The PSPC spectral index for both quasars is 2.3 ± 0. 1 (Fiore 
et al. 1994), steeper than any of the values in Table 3. For 
NAB 0205 4 024 the ASCA low-energy index is flatter than the 
PSPC by 0.4-0.8, and the ASCA best-fitting W H coincides with the 
Galactic value, suggesting that the spectrum continues to steepen 
below ASCA X-ray energies. 

To parametrize a curved spectrum we have also fitted the 0.6-10 
keV spectra of the two quasars with: a broken power-law model; a 
power law plus optically thin plasma emission (Raymond & Smith 
1977) model; a power law plus blackbody model; a power law plus 
Compton reflection model. The results are again shown in Table 3 
and confirm the presence of significant curvature in the spectra of 
these quasars. 

In PG 1244 4 026 the power law plus Raymond- Smith model 
gives a x 2 significantly better than the broken power-law model 
because it also fits the ~l-keV feature (see next subsection). 
However, the best-fitting energy index is still very steep and positive 
residuals are evident above 4 keV (Fig. 6). In NAB 0205 4 024 the 
power law plus Raymond- Smith model gives an acceptable \ °nly 
for low metal abundances, given the absence of significant features 
in the spectrum, similar to the ROSAT results of Fiore et al. (1994). 
Above 4 keV the best fit is similar to that in the broken power-law 
model. Fits with a power law plus blackbody model give acceptable 
X 2 in both cases. The power-law indices in this case are slightly 
steeper than in the broken power-law model. Inspection of the 
residuals shows again a slight excess of counts at high energy. In 
PG 1244 4 026 the power law plus blackbody model again gives a 
X 2 significantly better than the broken power-law model because it 
partly fits also the feature around — 1 keV (see next subsection). Fits 
with a Compton reflection model (plrefl in xspec) give x sig- 
nificantly higher than the previous models and push the parameter 
12/2 it to high and implausible values. 

We do not see any significant Line emission at the energies of 
the iron Ka lines. The 90 per cent upper limits to the equivalent 
width of a narrow line at 6.4 (and 6.7) keV, rest frame, in 


ASCA SIS - GIS 


PG1 244+026 



Figure 6. The SIS 4 GIS spectra of PG 1244 4 026 fitted with a power law + Raymond- Smith model. The lower panel shows die ratio between the data and the 
model. 
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0 

J PG 1244 + 026 and NAB 0205 + 024 are 400 and 230 eV (640, 
314 eV) respectively. 

oj The x 2 for the fits to the PG 1 244 -f 026 spectra are much higher 

| than those for the similar fits to the NAB 0205 + 024 spectra 

1 because of the presence of the 4 1-keV' feature in the former 

<£ quasar. We discuss this feature next. 

rH 

4J2 Low-energy features 

We have already noted that a ‘ 1 -keV’ feature seems to be present 
only in PG 1244 + 026, the quasar with the lower luminosity. Fig. 7 
shows the results of a simple power* law model fit to the 0.6-4 keV 
spectrum of PG 1244 4- 026: the 4 1-keV* feature is clearly visible. 

It has recently been realized that Residual Dark Distribution 
(RDD; the error in ASCA * s onboard correction for CCD dark 
current) can affect the low-energy ASCA SIS spectra. The symptom 
of the RDD problem is a sudden decrease in the SIS effective area 
towards lower energies from about 1 keV. We arc however confident 
that this has little effect on the 1-keV’ feature, for the following 
reasons: (a) SIS data of PG 1244 *f 026 and NAB 0205 + 024 
have similar RDD value but the 4 1-keV’ feature is visible in 
the spectrum of PG 1244 + 026 only. The 90 per cent upper 
limit to any unresolved Gaussian line emission at about 1 keV in 
NAB 0205 + 024 is of 20 eV. (b) The feature is present also in 
the GIS spectrum (c). The feature is present at the same level in 
SIS0 and SIS1 T while the RDD tends to be greater for SIS1 than for 
SIS0. 

To investigate the physical origin of the ‘1-keV’ feature and to 
quantify its strength, we performed a series of spectral fits to the SIS 
and GIS spectra of PG 1 244 -f 026 using only the 0.6-4 keV range, 
to avoid the complications of the hard tail (see previous section). 
The results are reported in Table 4, along with those of the other fits 
in this section. 


4.2.1 Emission-line model 

The fit with a power-law model plus a Gaussian line gives a small x 2 
(90.1 for 101 degrees of freedom). The equivalent width of the 1- 
keV feature, in the power law plus Gaussian line fit, is similar in the 
SIS and GIS detectors. In Table 4 we report the SIS determination of 
64 ± 15 eV. The Iinewidth is well constrained and cleanly resolved 
in the SIS spectrum to o = 0. 1 1 to ol keV. 

A fit with a power law plus a thermal plasma model (Raymond & 
Smith 1977) gives an acceptable x 2 (93. 1 for 100 d.o.f.). The best-fit 
temperature (1 keV) implies that the 4 1-keV emission is dominated 
by a blend of iron L and neon emission lines. 

An emission-line feature can be mimicked by fitting a spectrum 
with a strong absorption feature at slightly higher energies. We have 
then fitted the SIS and GIS spectra with models including absorp- 
tion structures. 


4.2.2 Absorption-edge Jits 

Fits with one absorption edge give x 2 significantly higher than the 
previous case (102.5 for 102 d.o.f.). The edge energy (1.17 keV, 
quasar frame) is consistent with that of Ne ex and/or Fe L, Fe xvi and 
Fe xvu. The best-fitting neutral Nu is significantly higher than the 
Galactic value. This is reasonable, since if there is highly ionized Ne 
and FeL absorption it is likely to have also highly ionized oxygen 
absorption at 0.74-0.87 keV. We then refitted the SIS and GIS 
spectra with a model including three absorption edges, at the 
energies of the most abundant ions in highly ionized gas with 
high Neix abundance, i.e. Ovn, Ovra and Ne ix-Fe xvi-Fe xvn 
(Nicastro et al. 1998), fixing the cold // H to the Galactic value. The 
results were not satisfactory. The depth of the oxygen edges is zero 
with small upper limits and x 2 is significantly higher than in the 
previous case: 1 14.7. Leaving free improves x 2 > but the oxygen 


ASCA SIS - GlS 


PG 1 244+026 



chcmael energy (keV) 


Figure 7. The SIS + GIS spectra of PG 1 244 + 026 fitted in the 0.6-4 keV range. The lower panel shows the residuals when the spectrum is fitted with a simple 
absorbed power law. 
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Table 4. PG 1244 + 026: spectral fits in the band 0.6-4 keV. 


Continuum models 


model 



T b 

A c 


X 2 (dof) 

PL 

1.7010.02 

1. 9+1.0 

- 


- 

131.8 (104) 

PL+BB 

1.50±0.11 

9.515.0 

0.1410.02 

- 

- 

109.5 (101) 

PL+Raym 

1.6110.03 

1. 9+0.3 
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- 
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emission line 



line E fa 

line EW<* 


X 2 (dof) 

PL+Gauss 
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1.9+0. 5 
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0.2610.08 

- 
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U or logT 

z abs. 

X 2 (dof) 

1-Phot. ion. 

1.92±0.10 

11.113.5 

22.5910.15 

2.0010.12 

0.2510.07 

93.2 (101) 

2-Phot. ion. 

1.94±0.10 

11.013.5 

21.5610.10 
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-0.3310.10 

95.7 (101) 

3-Coll. ion. 

1.90±0.10 

11.313.5 

22.6210.14 

6.9210.07 

0.2410.06 

90.7 (101) 

absorption notch models 


T b , notch £ b 

notch width** 

Cov. frac. 

X 2 (dof) 

PL+notch 

1.93±0.08 

9.313.3 

1.4410.05 

680170 

0.1410.04 

97.5 (100) 

PL+BB+notch 

1.48±0.20 

7.814.4 

7=0.1610.03 



97.1 (100) 




1.1610.03 

1417 

1FIXED 



J In 10 20 cm -2 . b In keV. c Metal abundances. d In eV. 

edge depths are still zero and the fit resembles completely the 
single-edge fit. Fixing to the Galactic value but leaving free the 
energies of two edges again produces a good fit (x 2 = 94.8, 
101 d.o.fi). The best -fitting energy of one edge is again 1.12 keV, 
but that of the other edge is < 0.64 keV (observer frame), close to 
the lower boundary of the observed range. So, if the cold absorption 
is fixed to the Galactic value, then there must be additional 
absorption edge(s) at energies lower than the observed range, 
corresponding to oxygen less ionized than O vi. We note, however, 
that this conclusion is weakened by the unknown contribution of the 
SIS RDD, which pushes low-energy events below the detection 
threshold. 

A 1 . 1 7-keV absorption feature can also be interpreted in terms of 
blueshifted oxygen absorption (Leighly et al. 1997). In this case, 
assuming that the absorption is mostly due to Ovra, the shift from 
the quasar frame would be equivalent to z — —0,38. A more 
complex continuum has little effects on the best-fit parameters of 
absorption edges. 

4. 2. 3 Ionized absorber model fits 

We fitted the data with a detailed ionized absorber model (not 
including resonant scattering absorption lines). We first generated a 
grid of photoionization equilibrium models using cloudy (Ferland 
1996), and fitted the spectrum interpolating by this grid, using the 
method of Fiore et al. (1993). To calculate the models we have 
assumed the observed spectral energy distribution (SED) (Fiore et 
al. 1995; Elvis et al. 1994). This is important, since the soft X-ray 
spectrum of this source strongly differs from that of ‘normal* 
Seyfert 1 galaxies, where warm absorbers are usually found 
(Reynolds 1997). A steep soft X-ray spectrum can completely 
ionize oxygen and neon but not iron, and so can produce edges in 
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Energy (keV) 


Figure 8. The thick line shows the best-fitting photoionized absorber model 
for PG 1 244 + 026 using the observed steep spectrum (note that the absorber 
has z = 0.25, while the quasar has z = 0.048), The thin line shows a 
photoionization model obtained using a much flatter SED (see text). Both 
models produce Fe L absorption features around 1 keV The flat SED also 
produces a deep Ovm edge, which is not observed. 

the 1-2 keV (Fe L) and 7-9 keV (Fe K) ranges, but not in the 
‘oxygen* 0.6-0.9 keV band. A fit with this model produces an 
acceptable x 2 (see Table 4, ionized absorber model 1). 

In Fig. 8 we show the best-fitting steep SED model (thick line) 
and a photoionization model obtained using a standard, much flatter 
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AGN SED (a power law of a = 1.2 from ultraviolet to X-rays, 
thin line), which, above 1 keV, gives a comparably good fit to the 
data. While iron in the flat SED model has an ionization structure 
similar to that of the model obtained using the right SED, oxygen is 
much less ionized: note the deep O vra edge present in the flat SED 
model. 

In this fit the redshift of the absorber is significantly higher than 
that of the quasar, because the main feature in the transmitted 
spectrum is the Fe xvm edge at 1.36 keV, while the deepest edge in 
the quasar spectrum is at 1 . 1 7 keV (quasar frame). However, a good 
fit can also be obtained for a different absorber redshift ( z = —033, 
Table 4, ionized absorber model 2). In this case the 1 J 7 keV feature 
is interpreted in terms of Ovn and O vm absorption. We cannot 
discriminate between these two solutions on statistical grounds. 

We have also tried fits with a collisional equilibrium model 
(Nicastro et al. 1998). The results were very similar to those 
obtained in the case of photoionization equilibrium (see Table 4, 
ionized absorber model 3). 

In all fits with detailed warm absorber models the column of cold 
gas is significantly higher than the Galactic value, similar to the 
values found using single edges to parametrize the absorber 
(Table 4). 

4.2.4 Absorption-line fits 

Fits with a single Gaussian absorption line do not give acceptable 
X 2 . Fits with two or more absorption Gaussian lines can produce x 2 
of 96 or smaller These models are indistinguishable for the SIS 
from models with a broad absorption notch, which we discuss in the 
following. Physical models including resonant absorption lines, as 
well as absorption edges from ionized plasma, will be discussed in a 
future paper (Nicastro et al , in preparation). 

4.2.5 Absorption-notch fits 

Fits with an absorption notch give x 2 higher than those with an 
emission line by Ax 2 ~ 7. While the power law plus notch lit 
is formally acceptable, the best-fitting value of the notch width is 
implausibly large (almost 1 keV), forced by the very low value of 
the covering fraction required by the fit. However, the notch best- 
fitting parameters are strongly dependent on the proper modelling 
of the continuum. For example, using a power law plus a blackbody 
to parametrize the continuum gives an acceptable fit fixing the notch 
covering fraction to 1, which in turn results in a much more 
reasonable value for the notch width of 14 ± 7 eV. The ASCA 
bandwidth is not wide enough and its spectral resolution is not good 
enough to constrain adequately both a complex continuum and the 
notch parameters. 

4.2.6 Comparison with the PSPC results 

PG 1244 + 026 was observed with the PSPC in 1991 December and 


the results of this observation have been reported by Fiore et al. 
(1994). The 0.6-2 keV flux level during the PSPC observation was 
about 50 per cent lower than the mean flux in the ASCA observation. 
The fit of an absorbed power law to the PSPC spectrum gives an 
acceptable x 2 (22.8 for 25 d.o.f.), kev = 2.3 ± 0.1 and 
= 2.9 ± 0.3 x 10 20 cm” 2 , slightly higher than the Galactic 
value. Most of the PSPC counts were detected below 0.3 keV, in 
the ‘carbon’ band. The quality of the spectrum between 0.3 and 
2 keV is not very high and emission or absorption features fainter 
than ^20 per cent cannot be excluded in this energy band. No 
evidence of spectral variability is present in the ROSATdata despite 
a factor of 2 flux variability. 

The PSPC data strongly constrain the level of any cold or warm 
absorption affecting the ‘carbon* band. Best-fitting models to the 
ASCA data including absorption features in the 1-2 keV band 
(Table 4) require a rather large absorption in addition to the Galactic 
one below 1 keV. Therefore, it is important to study whether the 
ASCA best-fitting models are consistent with the PSPC ones. Rather 
than performing joint fits to the ASCA and PSPC data, which are 
complicated by the large uncertainty in the relative PSPC/ASCA SIS 
calibration, and by the detailed shape of the continuum over the 
broad 0.1-4 keV band, we fitted the PSPC data with a power-law 
model including the emission or absorption features found in the 
previous section (see Table 4). The results are shown in Table 5, 
where we also report (in brackets) the 99 per cent parameter upper 
limits, or confidence intervals, when appropriate. 

We see that the presence of an emission line at 0.91 keV is not 
required by the PSPC spectrum, but an equivalent width of 65 eV is 
not excluded (10 per cent probability). However, the presence of a 
cold absorber of thickness 8.2 ± 3.0 x 10 20 cm -2 is inconsistent 
with the PSPC spectrum (probability < 1 per cent), while the 
presence of an edge at 0.62 keV with t = 0.51 ±0.14 is only 
marginally consistent with the PSPC result. 

5 DISCUSSION 
5.1 Continuum 

The ASCA observations of PG 1244 + 026 and NAB 0205 + 024 
have shown that the X-ray continuum of these two quasars flattens 
by 0.5-1 on passing from the 0.1-2 keV (PSPC) to the 2-10 keV 
band. Similar results were obtained by Brandt et al. (1997); 
by Comastri et al. (1998) and Leighly et al. (in preparation) on 
TONS 180, by Pounds et al. (1995) and Fiore et al. (in preparation) 
on REJ 1034 + 390, and by Leighly et al. (in preparation) on 
AKN 564. It appears that the X-ray spectrum of a sizeable 
number of steep PSPC and narrow Balmer line quasars has 
significant curvature, being flatter at higher energies. 

This could be a result of different components influencing the 
spectrum at different energies, as might happen in ‘normal’ broad- 
line Seyfert 1 galaxies and quasars, where a soft component is often 
present. The relative intensity of the two components would be 


Table 5. PG 1244 + 026: spectral fits to the PSPC. 


model 


a E 

line or edge E 6 

line EW C or t 

x 2 (dof) 

PL 

2.9±0.3 

2.3+0. 1 



22.8 125) 

PL+gauss 

2.9±0.3 

2.3±0.1 

0.91FIXED 

0+68 {+1001 

22.8 (23) 

PL+1 edge 

2.7±0.4(0.6) 

2.2±0.2 

1.12FIXED 

0.3^i 

20.9 (24) 

PL+2 edges 

2.7+0.4 

2. 2+0.2 

1.12F1XED 
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quite different from ‘normal’ quasars. Laor et al. ( 1 997) suggested a 
fainter hard component relative to the optical in the majority of low- 
redshift PG quasars (assuming a two-component model). However, 
Grupe (1996) found evidence for a stronger soft excess in a sample 
of soft X-ray selected Seyfert galaxies dominated by narrow -line 
objects. 

A large relative intensity of the soft component has important 
consequences on various competing models for the soft component. 
Disc reprocessing models (Matt, Fabian & Ross 1993; Fiore, Matt 
& Nicastro 1997) would require highly anisotropic emission to 
account for the discrepancy between the observed soft and hard 
fluxes. Optically thin free-free emission (e.g. Barvainis 1993) is 
ruled out by these observations because the best-fitting power-law 
slope is still too steep to fit the spectrum about 4 keV, because 
variability rules out optically thin plasma (Elvis et al. 1991) and 
because of the implausibly low metal abundances (< 1 per cent 
solar, see Table 3) required in NAB 0205 -f- 024 (see section 3.2. 1 ). 
The most likely origin for the steep component is Comptonized disc 
emission (e.g. Czerny & Elvis 1987; Fiore et al. 1995; Pounds et al. 
1995). 

The high-energy spectral index of the two quasars — 1.0 (see 
Table 3) is consistent with that of ‘normal’ Seyfert 1 galaxies (e.g. 
Nandra & Pounds 1994; Nandra et al. 1997). The enor on a H is 
however large and so no strong conclusion can be drawn on the 
origin of the hard emission. An answer to this question must await 
the large area and high energy sensitivity of AXAF, XMM and 
Spectrum X- gamma. 

5.2 Origin of the ~l-keV feature in PG 1244 4- 026 

The — 1-keV feature in PG 1244 4- 026 could be explained in terms 
of either a broad (o = 0.1 keV) emission line at 0.91 keV (0.95 keV 
quasar frame) of about 60 eV equivalent width or a r = 0.25 
absorption edge at 1.17 keV (or an absorption notch at 1.22 keV). 
These possibilities cannot be discriminated between on statistical 
grounds. 

An absorption interpretation requires additional low-energy 
absorption, either cold (with a column density higher than Galactic 
by ~~7 x 10 20 cm" 2 ; see Tables 4 and 5), in contrast with the 
PS PC results, or a warm absorber with a peculiar ion abundance 
distribution. 

The observed spectrum can be interpreted in terms of either an 
inflowing ( vie = 0.25) or an outflowing {vie = —0.33) absorber. In 
the first case the ion contributing most to the absorption is Fe xvm, 
and in the second case it is Ovm. The two cases cannot be 
discriminated between on statistical grounds. While an outflowing 
absorber has been suggested in several other cases (see e.g. Mathur 
et al. 1994), this would be the first case for an inflowing highly 
ionized absorber. 

A similar situation is found in IRAS 1 3224 — 3809 by Otani et al. 
(1996), and in AKN 564 by Brandt et al. (1994) and Leighly et al. 
(in preparation). Otani et al. (1996) and Leighly et al. (1997) 
interpret the features in the 1-2 keV band in terms of blueshifted 
absorption from relativistically (v = 0. 2-0.6 c) outflowing 
material. If the 1.17-keV feature in PG 1244 4 026 is due to a 
blueshifted Ovm photoelectric absorption, then the absorption seen 
below 0.64 keV may be due to C vi photoelectric absorption from 
the same gas. This distribution is far from an equilibrium distribu- 
tion (see e.g. Nicastro et al. 1998), not an impossible situation 
considering the large variability observed in this source. 

The absorption features seen between 1 and 2 keV may also be 
interpreted in terms of resonant lines (e.g. Leighly et al. 1997). If the 
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ion producing the absorption is oxygen Ovm (resonant line at 0.65 
keV), the best-fitting notch energy of 1.22 keV (quasar frame) 
implies a very high gas velocity: —0.56c. Fe xvm (E = 0.87 keV) or 
Fexvn (E = 0.81 keV) can also contribute to the absorption, 
because of their high oscillator strengths, 1.7 and 0.6 respectively 
(Kato 1976), and abundances. If the 1.22-keV absorption notch is 
caused by these ions, then the velocity of the outflowing gas will be 
smaller, —0.33c. 

In any case, in the blueshifted absorption scenario the gas is 
outflowing at velocities that are a sizeable fraction of c, reminiscent 
of blobs of gas in jets. It is interesting to note that similar absorption 
features have been observed in blazars (e.g. PKS 2155 — 304, 
Canizares & Kruper 1984; other BILacertae, Madejsky et al. 
1991; and 3C 273, Grandi et al. 1997), but usually below 1 keV. 
Somewhat surprisingly this implies less extreme conditions in these 
radio-loud objects than in our radio-quiet quasars. High-redshift 
radio-loud quasars may have similar jet-related absorption too 
(Elvis et al. 1997). 

An alternative interpretation of the 1 -keV feature is in terms of an 
emission line due to highly ionized oxygen, neon (Ne ix) and/or to 
iron L. There are two possible origins for this line: recombination in 
an optically thin thermal plasma, or ‘reflection’ in photoionized 
matter. 

5.2.2 Thermal plasma 

A thermal plasma is highly implausible on physical grounds (Elvis 
et al. 1991; Fiore et al. 1995). Emission measure is —2.5 x 10 65 
cm” 3 . For a spherical source (with radius R) and constant electron 
density n e , n t = 2.43 x 10 32 E” 3/2 cm” 3 . R< 10 14 cm from the 
observed X-ray variability, so n e < 1.5 x 10 12 cm” 3 . This implies 
an electron scattering optical depth r x 5: 30. With such values a 
thermal plasma is no longer optically thin. The situation is even 
worse if the matter is clumpy, as the density of each cloud must be 
greater. 

5.2.2 Photoionized matter 

The second possibility is that the ‘1-keV’ emitting marter is 
photoionized by the central nucleus. We can assume that the gas 
is not covering the source because there is no significant O vn and/or 
Ovm absorption. We therefore assume that we are not observing a 
‘warn absorber’, but rather a ‘warm reflector’. This could be either 
a warm absorber viewed from its side, in which case the matter 
would be optically thin to Thomson scattering; or the accretion disc, 
and the matter would be thick. Since both the ‘reflector’ and the 
primary emission (which provides most of the continuum) are 
observed, the optically thick case gives the highest values of the 
equivalent width (EW). However, even in the optically thick case, 
the expected EWs can barely account for the observed values (see 
below); hence we neglect the optically thin case altogether. 

The observed line (which is significantly broad) may be a blend 
of several lines (see for instance Zycki et al. 1994; Netzer 1996). 
The observed energy suggests the 0.92-keV Neix and 1.02-keV 
Nex recombination lines, the 0.87-keV Ovm recombination (to 
ground state) line and the Fe l (around 0.8 ke V) lines being the most 
important. None of these lines alone can account for the observed 
EW: for instance, the maximum value (i.e. for a face-on disc with an 
intervening ion fraction —0.6) for the Neix line is about 10 eV, 
while that of the oxygen recombination line is about 15 eV (note 
that an Ovm Ka recombination line at 0.65 keV with a similar EW 
should also be present; the 90 per cent upper limit on such a line is 
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30 eV). These values have been calculated using the formulae of 
Basco (1978), and assuming a reasonable ionization structure. (In 
the disc hypothesis a Fe Kar line at 6.5-6. 9 keV is also expected, but 
the upper limit of 300-400 eV does not exclude the presence of 
such a line.) Allowing for a possible factor of 2 neon and/or iron 
overabundance (an oxygen overabundance would decrease the Ne 
line while not increasing the O line) and/or anisotropy of the 
illuminating radiation, the observed equivalent width could be 
explained. (Note that in many Seyfert 1 galaxies the Fe Ka line is 
also stronger than expected, suggesting iron overabundance or 
anisotropic illumination.) 

In this scenario the ‘1-keV’ feature may arise from a highly 
ionized accretion disc. In the Matt et al. (1993) models, high 
ionization mainly arises because of high accretion rates (the 
ionization parameter depends on m 3 ). The detection of these 
emission lines in steep soft X-ray quasars would then be further 
evidence of high L/L&u. We note that, because recombination can 
occur only in highly ionized atoms, we would not expect features of 
this kind in ‘normal' quasars, as their disc should be much less 
ionized, as in fact observed. 

The strong dependence of the ionization parameter on m allows 
for large differences in the ionization structure against small 
differences in m and therefore that the k 1 -keV’ feature may not be 
ubiquitous in NLSyls. Indeed, a similar feature is not present in 
NAB 0205 + 024. A line with the same equivalent width as in 
PG 1244 + 026 (60 eV) would have been detected in the 
NAB 0205 + 024 SIS spectrum (the 90 per cent upper limit is 
only 20 eV). In the accretion disc scenario this would imply a 
different ionization state (higher or lower) of the matter or a 
significant metal undeTabundance in NAB 0205 4- 024. 

53 Variability 

The mean 0.5-10 keV luminosity measured by ASCA in 
PG 1244 + 026 and NAB 0205 + 024 differs by an order of 
magnitude: 7.3 x 10 43 and 6.4 x 10 44 erg s~ ! respectively (assum- 
ing isotropic emission, H 0 = 50.0 and <jo = 0.0). The optical 
(3000 A) monochromatic luminosities differ even more: 4x 10 43 
and 9.4 x 10 44 erg s -1 respectively. The variations seen in the NAB 
0205 + 024 light curve imply an efficiency in the conversion of 
matter into radiation greater than 1.6 per cent (e.g. Fabian 1984). 
Taken at face value, this excludes thermonuclear reactions as the 
origin of the observed X-ray luminosity for which the upper limit on 
the efficiency is in this case about 0.7 per cent 

The similar variability observed in the two quasars agrees with 
the Fiore et al. (1998) finding that the variability properties of (PG) 
quasars are correlated with the shape of the soft X-ray spectrum and 
the width of the Balmer lines, and so possibly then with the 
accretion rate, in the scheme of Pounds et al. (1995) and Laor 
et al. (1994, 1997). 

6 CONCLUSIONS 

ASCA observations of two steep soft X-ray quasars have shown that: 

(i) The X-ray continuum of the two quasars flattens by 
Aa = 0.5-1 going towards high energies. Similar results were 
obtained by other authors on some half-dozen steep soft X-ray 
quasars. 

(ii) PG 1244 + 026 shows a significant feature in the ‘1-keV’ 
region. Similar features were again reported in other steep soft X- 
ray quasars. The data are not good enough to discriminate between a 


broad emission fine centred at 0.95 keV (quasar frame) or an 
absorption edge at 1.17 keV, or an absorption notch at 1.22 keV. 
Line emission could be due to reflection from a highly ionized 
accretion disc, in line with the view that steep soft X-ray quasars are 
emitting close to the Eddington luminosity. Photoelectric edge 
absorption or resonant line absorption could be produced by gas 
outflowing at a large velocity (0.3-0.6 c). In these absorption 
models significant cold (i.e. oxygen less ionized than Ovi) absorp- 
tion in excess of the Galactic is required. This would imply an 
increase by a factor of 2-3 of the cold column with respect to a 
previous PSPC observation or a peculiar ionization structure. In 
neither the emission nor absorption cases is the SIS resolution good 
enough to identify unambiguously the ions responsible for the 
feature. The high resolution and high throughput of the low- 
energy gratings and spectrometers of AXAF and XMM are clearly 
needed to shed fight on this puzzling case. 

(iii) The two quasars show similar variability properties (flux 
variations up to a factor of 2 in 10 ks) despite a factor of 10 
difference in the X-ray observed luminosity. This agrees with the 
Fiore et al. (1998) finding that the variability properties of radio- 
quiet quasars are correlated with the shape of the X-ray spectrum, 
the width of the Balmer lines and so possibly with the accretion rate. 
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We present timing and spectral analyses of the ASCA observation of the Seyfert 1 galaxy NGC 985. The 
0.6-10 keV spectrum of this source is complex: large residuals are evident below 1 keV when fitting the spectrum 
with a power-law model. Fitting a warm absorber model to the 0.6-2. 5 keV spectrum gives a = 1.12 10.04, 
LogA r ^ M = 21.97 1 0.08 and LogU= 0.061 0.09, but the residuals continue to show a deficit of counts between 0.9 
and 1 keV. Adding an absorption line improves the fit, and the energy of the line is consistent with that of Ka 
NeIX-X resonant absorption lines. Hence, we confirm the presence of an ionized absorber along the line of sight to 
this source and interpret the further 1 keV spectral feature as the first detection of a strong resonant absorption 
line associated with this system. The extrapolation of this model above 2.5 keV produces large positive residuals 
above 3-4 keV. Fitting the data with a broken power law plus warm absorber model gives an acceptable y 2 and 
Aq ~ 0.5. A narrow iron line at 6.4 keV (quasar frame) of equivalent width 1381 eV is also present in the 
ASCA data. 


1. Introduction 

Absorption features from ionized gets are com- 
mon in the X-ray spectra of type 1 AGN ([!]). 
Brandt et al. (1994, [2]) found a deep ab- 
sorption edge in the ROSAT-PSPC spectrum of 
NGC 985, at an energy consistent with OVII- 
OVIII K edges. 

Absorption K and L edges are not the only ab- 
sorption features expected to be imprinted on X- 
ray spectra of warm absorbed sources, if the gas 
is not isotropically distributed around the ioniz- 
ing source, resonant absorption lines should be 
present as well ([3]). Until now, none of the pre- 
dicted resonant absorption lines have been clearly 
identified in the X-ray spectra of AGN known to 
host warm absorbers. 

We present here the analysis of an ASCA obser- 
vation of NGC 985. We confirm the presence of a 
warm absorber in this source. We also detect an 
additional absorption feature at 0.97 keV, which 
could be explained in terms of resonant scattering 


absorption from the warm gets. 

2. Data Reduction and Analysis 

NGC 985 was observed by ASCA on August 
12, 1996 with a total exposure of ~ 40 ksec. SIS 
and GIS count rates were 0.25 counts s^ 1 and 
0.18 counts s^ 1 respectively. 

We applied standard screening procedures (fol- 
lowing Reynolds 1997, [1]). We extracted source 
counts from regions of 4’ and 6’ radius for SIS 
and GIS respectively. Both local and blank-sky 
field backgrounds were tested. We did not find 
significant differences between the two cases. In 
the following we present the results obtained us- 
ing the background from the blank fields. 

2.1. Timing Analysis 

The source intensity changes by a factor ~ 1.4 
during the observation. Fig. 1 shows the SIS 
and GIS light curves. We grouped the data in 
bins of 500 s. Changes are more marked in the 
SIS light curve, suggesting spectral variability as- 
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Time (in 10 4 s) 

Figure 1. SIS and GIS light curves of the ASCA 
observation of NGC 985. 


sociated with the soft component. However the 
quality of the data and the small dynamical range 
of variability does not allow a detailed study of 
the spectral variability. 

2.2. Spectral Analysis 

We fitted the SIS spectra in the 0.6-8 keV en- 
ergy band and the GIS spectra in the 1-10 keV 
band. 

2.2.1. Soft (0.6-2. 5 keV) X-ray Spectrum 

Table 1 shows the results of the fits to the 0.6- 
2.5 keV SIS and GIS spectra. In Fig. 2 we show 
the ratio between data and best fit models. We 
first fitted the spectra with a simple power law. 
reduced at low energy by cold absorption (Model 
A). This fit is unacceptable and the absorption is 
not consistent with the Galactic value along the 
line of sight (3.0xl0 20 cm -2 ). Large positive and 
negative residuals are evident below ~ 1 keV in 
the upper panel of Fig. 2, suggesting the presence 
of either a soft excess, a warm absorber compo- 
nent. or both. A warm absorber component asso- 
ciated with NGC 985 has already been proposed 
by Brandt et al. (1994, [2]) in their analysis of 
ROSAT data of this source. We therefore added a 
photoionization model (built with CLOUDY, ver. 

90.01, [4]) to Model A, and refitted the data fixing 
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Figure 2. Residuals after fitting the 0.6-2. 5 keV 
data with models A (upper panel), B (second 
panel) and C (third panel). The lower panel 
shows the residuals after the extrapolation of the 
0.6-2. 5 keV best fit Model C down to 0.45 keV. 


Nh to the Galactic value (Model B). Power-law 
index, U and N^ A were left free to vary. Residu- 
als to the best fit model are shown in the second 
panel of Fig. 2: a deficit of counts is still evident 
between 0.9 and 1 keV. 

We then added a negative Gaussian to the previ- 
ous model, and refitted the data (Model C). The 
addition of this component provides a highly sig- 
nificant improvement in y 2 : F(2,40) = 11.6, cor- 
responding to a probability of < 1%. Residuals 
to this model are shown in the third panel of Fig. 
2. The best fit continuum and warm absorber pa- 
rameters are quite similar to those obtained fit- 
ting the data without the absorption line. The 
best fit line energy is 0. 97^^05^ consistent with 
a blend of NeIX-X Ka resonant absorption lines 
(at the source redshift of 0.043). The equivalent 
width is 18.7 eV. 

2.2.2. 0.6-10 keV X-ray Spectrum 

Fig. 3 shows the residuals from the extrapola- 
tion of the best fit Model C up to 10 keV. A clear 
flattening of the spectrum is present above 3 keV. 
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Figure 3. Residuals after the extrapolation of the 
best fit Model C up to 10 keV. A clear flattening 
of the spectrum is present above 2-3 keV. 


We then repeated the fit substituting in Model C 
the simple power law with a broken power law, 
fixing the low energy spectral index to the best 
fit value of Model C, and allowing the break en- 
ergy and the hard spectral index to vary (Model 
D). Results are presented in Table 2. 

Finally we added a narrow line at 6.4 keV (quasar 
frame), to Model D, and refitted the data (Model 
E: table 2). The addition of the iron line to Model 
D is significant: F(l,108)=7, corresponding to 
P~l%. The equivalent width is 138^ i 0 eV. 

3. Physical Interpretation of the 1 keV 
spectral feature 

A number of possible interpretations could be 
invoked to explain the 1 keV absorption fea- 
ture present in the ASCA spectrum of NGC 985. 


These include non-equilibrium effects, which 
would produce emerging spectra quite different 
from those expected at equilibrium ([5]), multi- 
zone geometrical configurations, in which a num- 
ber of clouds of gas with different column densi- 
ties and ionization states would obscure the line 
of sight and imprint different absorption features 
in the emerging spectrum ([6]), or continuum ef- 
fects. Another possible explanation is provided 
by resonance absorption lines in the context of 
a simple single-zone, equilibrium photoionizaion 
model ([7]). In the following we discuss this lat- 
ter possibility. A complete presentation of all the 
above possibilities is deferred to a forthcoming 
paper. 

Resonant absorption lines are expected in the 
X-ray spectra emerging from ionized gas. Many 
of these atomic transitions have negligible oscil- 
lator strengths, and/or are produced by elements 
with low cosmic abundance. In a few cases (Kq 
resonant lines of C, O, Ne and Fe, and L iron 
lines), the lines are particularly intense (equiva- 
lent widths > 20 eV for turbulence velocity >100 
km/s) and have energies between 0.1 and 10 keV. 
We have built a simple code to calculate the 
equivalent width (and so the growth curves) and 
the profile of the main resonant absorption lines 
of C, O, Ne and Fe. The intensities of these lines 
depend on the geometrical configuration and dis- 
tribution of the gas around the central source. 
We run our program using the best fit Nff and U 
values of Model C, and assuming a simple geomet- 
rical configuration: a single cloud obscuring our 
line of sight. Fig. 4 shows the spectra (rest frame) 
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Table 2 

Spectral fits to the 0.6-10 keV SIS+GIS spectra 

Model qh “norm. ( ’EW em Xr(dof) 

D 2.7570.25 0.5170.06 3.3570.05 - 0.95(108) 

E 2.6670.25 0.5670.05 3.3570.05 138+^o 090(107) 

" In 10 -3 ph cm -2 s _1 keV -1 , at 1 keV. 4 5 In eV. 

note: We list only the parameters left free to vary in the fit. Other parameters 
are fixed to the best fit value of Model C. The emission line energy is fixed to 
6.4 keV (source frame). 


transmitted from such a cloud of ionized gas (the 
edge-like absorption features are those predicted 
by CLOUDY), for turbulence velocities of 100 
km/s (left panel) and 1000 km/s (right panel) 
respectively. The dashed line at 0.6 ke\ shows 
the lower end of the SIS band. OVII and OVIII 
absorption K edges are clearly visible in the spec- 
tra, along with the strong Ka resonant absorp- 
tion lines of NeIX and NeX. The total predicted 
equivalent width of the NeIX-X lines ranges from 
5 eV to 20 eV, for turbulence velocities between 
100 and 1000 km/s. OVII and OVIII Ka reso- 
nant absorption lines are also present in the spec- 
tra of Fig. 4, with total equivalent widths of 5.5 
eV (left panel) and 30 eV (right panel). They 
lie at the lower end of the considered SIS energy 
range. However residuals after the extrapolation 
of the 0.6-2. 5 keV best fit Model C down to 0.45 
keV show an absorption-like feature (Fig. 2, lower 
panel). The uncertain SIS calibration below 0.6 
keV does not allow us to draw a definitive con- 
clusion. 

4. Conclusion 

We present timing and spectral analyses of the 
ASCA data of the Seyfert 1 galaxy XGC 985, and 
find the following main results: (a) we confirm 
the presence of a warm absorber along the line 
of sight of NGC 985; (b) we detect an additional 
absorption feature at 0.97 keV (rest frame); (c) 
we interpret this feature as due to a blend of NeIX 
and NeX Ko resonant absorption lines from the 
ionized gas; (d) the spectrum of NGC 985 flattens 
above 3 keV by Aa = 0.5; (e) a narrow emission 
line at 6.4 keV (quasar frame) is also seen in the 
data, with an equivalent width of 138 eV . 



Energy (in keV) 

Figure 4. Spectra (in the observer frame) trans- 
mitted from a cloud of ionized gas with N H and 
U from the best fit values of Model C, for turbu- 
lence velocities of 100 km/s (left panel) and 1000 
km/s (right panel) respectively. 
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ABSTRACT 

We present photoionization models accounting for both photoelectric and resonant absorption. Reso- 
nance absorption lines from C, O, Ne, Mg, Si, S, and Fe between 0.1 and 10 keV are treated. In particu- 
lar we consider the complex of almost 60 strong Fe L absorption lines around 1 keV. We calculate 
profiles, intensities, and equivalent widths of each line, considering both Doppler and natural broadening 
mechanisms. Doppler broadening includes a term accounting for turbulence of the gas along the line of 
sight. We computed spectra transmitted by gas illuminated by drastically different ionizing continua and 
compared them to spectra observed in flat X-ray spectrum, broad optical emission-line type 1 active 
galactic nuclei (AGNs), and steep X-ray spectrum, narrow optical emission-line type 1 AGNs. We show 
that the ~ 1 keV absorption feature observed in moderate-resolution X-ray spectra of several narrow-line 
Seyfert 1 galaxies can be explained by photoionization models, taking into account resonance absorp- 
tion, without requiring relativistic outflowing velocities of the gas, if the physical properties of these 
absorbers are close to those found in flat X-ray spectrum Seyfert 1 galaxies. We finally present simula- 
tions of the spectra emerging from gas illuminated by both steep and flat ionizing continua, as seen by 
the AXAF high-energy transmission gratings and the baseline Constellation-X calorimeter. We discuss 
briefly the relevant physics that can be investigated with these instruments. 

Subject headings: galaxies: active — galaxies: Seyfert — line: formation — radiative transfer — 

X-rays: galaxies 


1. INTRODUCTION 

The ROSAT (Trumper 1983) PSPC (Pfefferman et al. 
1987) has found a large spread in the spectral indices of type 
1 active galactic nuclei (AGNs): 0.5 < a E <3 [/(£) = E~ XE ; 
Laor et al. 1997; Walter & Fink 1993; Fiore et al. 1994], 
This remarkably broad soft X-ray spectral index distribu- 
tion appears not to be random, but rather strongly corre- 
lated with other important optical and X-ray properties 
of these objects. Objects with narrower H/3 
(1000 < FWHM < 2000 km s _1 ) and smaller [O m]/ 
H/?(< 1) ratio are systematically steeper in soft X-rays 
(a £ > 1.5) than AGNs with broad (FWHM > 2000 km s _ *) 
permitted lines and greater [O m]/H/J ratio (1 < [O m]/ 
H/?<3, apspc < 1.5). In particular most of the so-called 
narrow-line Seyfert 1 galaxies (NLSyls) show a steep 0.2-2 
keV continuum (Boiler, Brandt, & Fink 1996). In the 2-10 
keV energy range narrow optical emission-line quasars 
have still a somewhat steeper spectrum compared to broad 
optical emission-line quasars. The spread of the 2-10 keV 
spectral index is, however, smaller than that of the 0.1-2 
keV PSPC one (Brandt, Mathur, & Elvis 1997). It has been 
suggested that these properties are correlated with the acc- 
retion rate (Pounds, Done, & Osborne 1995; Laor et al. 
1997; Fiore et al. 1998a, 1998b). 

Ionized gas, modifying the AGN X-ray spectrum, has 
been detected in about half of the bright, X-ray flat, broad 
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optical emission-line type 1 AGNs studied by ASCA 
(Reynolds 1997; George et al. 1998). The main features 
imprinted by the ionized gas on the X-ray spectrum are the 
O vu and O vin K absorption edges at 0.74 and 0.87 keV. 
As already pointed out by George et al. (1998), the simi- 
larity of the parameters of these features among different 
sources in both the Reynolds (1997) and the George et al. 
(1998) samples implies that “warm ” absorbers in AGNs all 
have roughly similar ionization structure, despite the broad 
range of ionizing luminosity (i.e., about 3 orders of magni- 
tude, from ~10 41 to 10 44 ergs s" 1 ). In terms of the ioniza- 
tion parameter U, the presence in the gas of He-like and 
H-like ions of C, O, and Ne as dominant ionic species 
implies a narrow range of log U = 0.5-1. 5 (Nicastro et al. 
1999). This is consistent with the location of the warm gas at 
radii scaling with the square root of the luminosity. The 
column densities of the Reynolds sample are distributed in 
the range log 7V H = 21—23, but the standard deviation of the 
distribution is rather small, <r log Nh = 0.3, again indicating a 
relative uniformity in the physical properties of this AGN 
component. Several type 1 AGNs showing ionized absorp- 
tion features in X-ray also show high-ionization UV 
absorption lines (namely, O vi, N v, C iv, C in]), superposed 
to the broad emission lines and systematically blueshifted 
with respect to them (Mathur et al. 1994; Mathur, Elvis, & 
Wilkes 1995; Mathur, Wilkes, & Elvis 1998). The widths of 
these absorption lines are typically much broader than the 
associated thermal broadening, suggesting turbulence 
velocities of the absorber. The UV and the X-ray absorbers 
have ionization states compatible with each other, which 
supports the identification of these two components 
(Mathur et al. 1994, 1995, 1998; Mathur, Wilkes, & Aldcroft 
1997). In one of these cases (the Seyfert 1 galaxy NGC 3516) 
the total column density implied by some strong UV 
absorption lines is larger than that required by the X-ray 
absorption features (Kriss et al. 1996a, 1996b). However it is 
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quite possible that more than one kinematic component 
contributes to the total intensity and equivalent width of 
the UV lines and that only one of those is physically associ- 
ated with the X-ray absorber (Mathur et al. 1997). The UV 
absorption lines may be produced in a turbulent outflowing 
wind intercepting the line of sight (possibly driven out by 
radiative acceleration). The outflowing velocities of the UV 
warm gas span the range v ouX ~ 100-2000 km s 1 (Mathur 
et al. 1994, 1995, 1997), and a similar range has been esti- 
mated for the dispersion velocities a v along the line of sight. 

On the other hand, deep absorption features at the ener- 
gies of the K edges of highly ionized oxygen have never been 
observed in the 0.1-2 keV spectra of soft X-ray steep spec- 
trum, narrow optical emission-line type 1 AGNs. Detec- 
tions of deep oxygen edges have been reported for two 
sources previously classified as NLSyl (IRAS 17020 + 4544, 
Leighly et al. 1997; Mrk 766, Matt et al. 1999). During the 
ASCA (0.5-10 keV) and BeppoSAX (0.1-10 keV) obser- 
vations the two sources have a rather fiat X-ray continuum 
(a £ ~ 1) and therefore cannot be considered X-ray steep 
AGNs. Large variations of the X-ray spectral shape have 
already been observed in the past (e.g., 1H 0419-577, 
Guainazzi et al. 1999, Turner et al. 1999; NGC 4051, 
Guainazzi et al. 1998; RE J2248 — 511, Puchnarewicz et al. 
1995). It is quite possible that a population of (perhaps low 
mass) AGNs does exist, switching between soft, high states 
and hard, low states (depending possibly on the accretion 
rate), and therefore classified in one way or in the other 
depending on the observed state (in analogy to what is seen 
in Galactic black hole candidates). 

Despite the lack of oxygen edges, steep soft X-ray spec- 
trum, narrow optical emission-line, type 1 AGNs are 
usually far from being featureless (Brandt 1995; Hayashida 
1997). ASCA has found large residuals between 1 and 2 keV, 
after subtracting the appropriate continuum, in at least four 
cases (Leighly et al. 1997; Fiore et al. 1998b). The interpre- 
tation of these features is not unique, given the relatively 
modest spectral resolution of the ASCA SIS and the com- 
plexity of the continua. Leighly et al. (1997) interpret the 
features found in the three NLSyls studied in terms of O vn 
and O vin resonance absorption lines blushifted by ~0.5c. 
Fiore et al. (1998b) interpret the feature found in PG 
1244 + 026 in terms of either an emission line at 0.92 keV or 
Fe xviii and Fe xvii resonant absorption blueshifted by 
~0.3c. 

Photoionization models used to fit and interpret AGN 
spectra usually deal with photoelectric absorption only. 
Matt (1994) and Krolik & Kriss (1995) were the first to 
point out the importance of X-ray resonant absorption lines 
in ionized absorbers. Here we calculate in detail photoion- 
ization models, including all the relevant resonant absorp- 
tion lines, and study the dependence of the transmitted 
spectrum on the exact shape of the ionizing continuum. We 
find that the phenomenologies observed in both broad 
optical emission-line type 1 AGNs and NLSyl galaxies can 
be explained in terms of their quite different ionizing contin- 
ua, provided that, for a given luminosity, the radial distance 
of the gas clouds from the ionizing source (in units of gravi- 
tational radii) and the gas density are similar in both classes 
of objects. 

Netzer (1993, 1996) discussed the relevance of recombi- 
nation and fluorescent emission lines, provided that the 
covering factor of the ionized matter is significant, whereas 
Krolik & Kriss (1995) and Matt, Brandt & Fabian (1996) 


pointed out that resonant scattering lines can be even more 
intense. Here, for the sake of simplicity, we do not consider 
any line emission, and therefore our model is strictly valid 
only if the covering factor is negligible (see next section). 

In § 2 we present our photoelectric plus resonant absorp- 
tion model and discuss the relevant physics. Section 2.3 
compares the transmitted spectra computed with our 
models for two different ionizing continua with the spectra 
of a broad optical emission-line type 1 AGNs and a NLSyl 
galaxy actually observed by ASCA (A E — 100 eV, effective 
area ~ 500 cm 2 ). Section 3 presents a set of simulations with 
instruments of much higher spectral resolution (AXAF 
high-energy transmission gratings [HETG] A E = 1.5 eV, 
collecting area - 100 cm 2 , AXAF Proposer's Guide version 
1.0, 1997 6 ; and Constellation-X baseline calorimeter, 
A£ = 3 eV, collecting area ~ 10,000 cm 2 . The High 
Throughput X-Ray Spectroscopy [HTXS] Mission, 1997 7 ). 
We finally present our conclusions in § 4. 

2. MODELS 

We calculated the curves of growth of all the 268 permit- 
ted resonant absorption lines of C, O, Ne, Mg, Si. S, and Fe 
in the 0.1-10 keV energy range with oscillator strengths 
greater than 0.05. We used atomic data from the complete 
list of permitted LS-type (LS coupling transitions) resonant 
lines of all ionization states of 18 elements tabulated in 
Verner, Verner, & Ferland (1996). Intercombination and 
forbidden lines are not included in that compilation. We 
carried out the calculations considering the exact Voigt 
profile for the lines, i.e., considering both Doppler and 
natural broadening. As suggested by the width of the 
absorption lines detected in the UV spectra of the UV/X- 
ray ionized absorbers (Mathur et al. 1994, 1995, 1997, 1998), 
in addition to the thermal Doppler broadening we consider 
here a term accounting for gas turbulence and refer to this 
term as the “ turbulence velocity ” a v along the line of sight. 

We calculated spectra emerging from photoionized gas 
for a number of values of N H , U, and a v . Relative ionic 
abundances and edge-like absorption features for given U 
and N h and for a particular shape of the nuclear ionizing 
continuum are those calculated by CLOUDY (Ferland 
1997, version 90.04), whereas the intensity and profile of the 
resonance absorption lines are consistently computed here 
for a given value of <r t .. The result is an emerging spectrum 
including both absorption edges and lines. 

We considered only line absorption, not emission. We 
used CLOUDY (version 90.04, Ferland 1997) to build total 
(transmitted plus emitted) spectra and verified that, for the 
range of total column densities and ionization parameters 
considered here, this assumption corresponds (for a spher- 
ical configuration) to covering factors f c < 0.1 (in good 
agreement with the Netzer 1993 results; see his Fig. 3). In 
an ideal optically thin case in which (1) the ionized gas 
isotropically and uniformly covers the central source ( f c = 
1), (2) any systematic (inflowing or outflowing) motion of 
the gas is absent, and (3) any line destruction mechanism 
(autoionization, Compton scattering, photoabsorption) is 
inefficient, then scattering of the resonant emission-line 
photons in the line of sight would exactly compensate for 
the line absorption along the line of sight, and, according to 

6 http ://asc. harvard.edu/udocs/docs/PG_frame.htnil. 

7 http ://constellation.gsfc.nasa.gov/documentation.html. 
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Netzer (1996), only emission features (mainly recombi- 
nation lines) would be present in the total spectra. Any 
deviation from the above ideal condition would produce net 
absorption lines, which can therefore diagnose the physical 
and geometrical conditions of the absorbing material. Pres- 
ence of outflows instead would imply that both emission 
and absorption lines are present in the emergent spectra, 
with the absorption lines blueshifted with respect to the 
emission ones. In this paper, however, for the sake of simpli- 
city we assume that a cloud of gas with negligible covering 
factor (f c < 0.1) is completely obscuring the line of sight, 
and therefore that the contribution from emission lines to 
the overall spectrum is negligible. In comparing our model 
with real spectra it must be kept in mind that absorption 
lines may be partly compensated by line scattering and 
recombination lines for significant covering factors or that 
emission lines may appear redward (blueward) of the 
absorption lines if the matter is outflowing (inflowing). 

In the next section we present the details of the computa- 
tion and show some examples of both growth curves and 
emerging spectra. 


2.1. Computation 

The equivalent width of a resonant absorption line is 


00 

EW = j (1 -e~ x ')dv y 


( 1 ) 


where t v is the dimensionless specific optical depth of the 
considered transition 


C L Tie 2 

t v = dsa v = n t L — . (2) 

Jo wc 

In the above equation a v is the absorption coefficient at the 
frequency v in cm -1 , n x is the number density of ions of the 
given element that populates the lower level in the unit 
volume (in cm ' 3 ), L is the linear dimension of the cloud of 
gas along the line of sight (in centimeters), / lu is the oscillator 
strength of the electron transition from the lower to the 
upper level, m is the electron mass in grams, and 0(v) is the 
normalized Voigt profile. With these units the equivalent 
width in equation (1) is in hertz. The Voigt profile can be 
written as (e.g., Rybicky & Lightman 1979) 


<P(v) = 


V / ^ Av D 

where H(a , a) is the Voigt function: 


H(a, u ) , 


a 

H(a , u) = - - 

* J-oo * 


e yl dy 
+ (u- y) 


2 * 


(3) 

(4) 


In equation (4) a is essentially the ratio between the natural 
and the Doppler line widths T and Av^, and u = (v 
— v 0 )/Av d . For small values of a the center of the line is 
dominated by the Doppler profile, whereas the wings are 
dominated by the Lorentzian profile. As is well known, this 
influences the growth curves, making them increase with JV H 
even over the core saturation. 

We consider two different contributions to the Doppler 
profile: (1) the thermal motion of the ions and (2) the disper- 
sion velocity <r v of the gas along the line of sight : 


A v D 


12 

c 




(5) 


Let us estimate the relative contribution of the two terms in 
equation (5): 


..therm 


J2kf 


~ 1.3 


m: 


05 A' 05 


m<r v 


K)i 


( 6 ) 


In equation (6) (T) 6 is the equilibrium temperature of the 
gas in units of 10 6 K, (Oxoo * s dispersion velocity in 
units of 100 km s -1 , and A is the atomic weight of the 
considered element. For a highly ionized gas in photoion- 
ization equilibrium, the temperature spans the range 
(T) 6 ~ 0.1-1. In such a gas the elements lighter than O-Ne 
are usually fully stripped (the threshold element over which 
bound ions are still present depending strongly on the pho- 
toionizing continuum shape), and so the residual opacity 
from both bound-free and bound-bound transitions is 
mainly due to the heavy elements (from O to Fe). Assuming 
therefore (7% = 0.5, (aJxoo = 1, and ^4 = 16 (the atomic 
weight of the lightest non-fully stripped element), we get, 
from equation (6), v thcrm /<r v ^ 0.2. For the iron ions (A = 56) 
we would find v thcrm /o v ~ 0.1. Hence, for dispersion veloci- 
ties of ~100 km s“\ the contribution of the thermal 
motion of the ions to the broadening of the absorption line 
is only — 10%-20% of the total, and it decreases linearly 
with increasing <j v . Therefore, the second term in equation 
(5) is the dominant one for (<r v ) 1 00 > 1. 


2.2. Growth Curves 

Figures 1 and 2 show the growth curves of the Ka and KjS 
lines of H-like and He-like C, O, Ne, and Fe, for (1) o v = 0 
(Doppler broadening entirely due to the thermal motion of 
the ions), and (2) o v = 1000 km s“ \ respectively. The equiv- 
alent widths of the two Ka lines of H-like ions have been 
summed to give a single equivalent width for the doublet. 

The ionization state of the absorbing gas is that of a 
typical flat X-ray spectrum, broad optical emission-line type 
1 AGN “warm absorber” (log U = 1) for C, O, and Ne, 
whereas a much higher ionization case (log U = 3) is 
adopted for the Fe xxv and Fe xxvi lines. In both cases the 
ionizing continuum used is the AGN spectral energy dis- 
tribution (SED), tabulated in CLOUDY (Ferland 1997, 
version 90.04). However, we stress that the growth curves 
depend only on the relative abundance of the ion under 
investigation and not on the particular combination of U 
and shape of the ionizing continuum producing this abun- 
dance. 

At low column densities (<J10 21 -10 23 cm 2 , the exact 
value depending on the considered line and dispersion 
velocity), the equivalent width of each line increases linearly 
with iV H . For higher values of N H the optical depth of the 
line core saturates and the Lorentzian profile begins to 
drive the further slower growth of the equivalent width with 
increasing columns. The linear branch of the growth curves 
(the low column density limit) is described by equation (1) in 
the limit r v < 1 : 


ne 2 

EW ^ dvz v = n,L — — f lu = r tot . (7) 

Jo mc 

Let us define A x , £ xt , and N H as the abundance of the 
element X relative to the H abundance, the relative abun- 
dance of the ion i of the element X, and the equivalent 
hydrogen column density along the line of sight, respec- 
tively. With these positions and introducing an appropriate 
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Fig. 1. — Growth curves of the Ka and K/S lines of H-like and He-like C, O, Ne, and Fe, for a v — 0. The ionization state 
typical warm absorber (log U = 1) for C, O, and Ne and a much higher ionization case (log U = 3) for Fe xxv and Fe xxvi. 


of the absorbing gas is that 


of a 


numerical format we have 

EW » 2.2(£ X i) 0 . 4 (Ax)_ 5 (N) 2 2(/iu)o.5 eV . (8) 

When the core of the line becomes opaque to the radi- 
ation, equation (8) is no longer valid, and the computation 
of the equivalent width must be carried out using equation 
(1) with the exact shape of the profile function. A compari- 
son between the curves in Figure 1 and those in Figure 2 
allows evaluation of the effects of the line profile upon the 
resulting equivalent width. When the Doppler broadening 
is dominated by the thermal motion of the ions (o v <t 100 
km s -1 , Fig. 1), the dynamical range of variation of the 
equivalent width of each line with increasing (V H is much 
smaller than in the case in which the Doppler broadening is 
dominated by a high dispersion velocity of the gas along the 
line of sight (a v > 100 km s~ x ). This is because when the line 
is strongly broadened the core of the line needs a larger 
column of gas to become optically thick, and so the linear 
branch of the growth curve spans over a larger range of N H 
before saturating. From equations (2), (3), and (5), using the 


same numerical factor as in equation (8), and setting {S at = 
ajc, we get 

to — 2-9(/ lu ) 05 (Cxi)o.4(4x )- 5 4 05 (N) 22 

x [(T) 6 + A( 2.3 x 10 3 j8J 2 ]-°' 5 £o 1 • (9) 

In the above equation the two terms within the square 
brackets describe the contribution to the core optical depth 
by the two considered Doppler mechanisms. From equa- 
tion (6) we know that dispersion velocities of the gas along 
the line of sight greater than ~ 100 km s _1 make the turbu- 
lence broadening mechanism very efficient compared to the 
thermal one and so reduce the core optical depth of each 
line. Equation (9) then allows definition and evaluation of 
the core optical depth in the two extreme cases: <x v = 0 
(thermal motion only, T' 0 herm ) and <r v S> ( kT/M ) 1/2 . We get 

~ 2.9(/ lu ) 0 . 5 (^ xi )o.4(4x)-5 A°- 5 (N) 22 (T) 6 -°- 5 £o 1 . 

(10) 

c b - 1.24(/Jo.s(^)o.4Ux)-5(N) 22 (^):^£o 1 • (11) 



112 


NICASTRO, FIORE, & MATT 


Vol. 517 


a; 

C 


-C 

TJ 


£ 

£ 

33 

> 

'5 

c r 
fcj 


> 

0) 

c 


cc 

> 


U 


21 22 23 24 21 22 23 24 



21 22 23 24 21 22 23 24 


LogN H 


L°gN H 


Fig. 2— Same as Fig. 1, but for a v = 1000 kms" 1 . 


The ratio between these two quantities does not depend on 
the particular considered transition, the distribution of the 
ionic abundances in the gas, and its column density, and is 
therefore a good estimator of the relative contribution of 
the two Doppler broadening mechanisms to the core 
optical depth : 

^therm 

= 2.3^° 5 (T)g °- 5 (^J_3 . (12) 

T o 

Equation (12) tells us that a given column of nonturbulent 
highly ionized gas in photoionization equilibrium can be 
more than a factor of 10 thicker to the resonant absorption 
process (at the energy of the transition) than an identical 
column of gas undergoing strong turbulence \o v > 300 km 

s-\ix.,0?J_ 3 >l]. 

2.3. Emerging X-Ray Spectra 

The relative importance of the spectral features in spectra 
emerging from photoionized gas obscuring the line of sight 
strongly depends on the spectral energy distribution of the 
ionizing continuum in the UV to X-ray energy range. We 


therefore treat separately the two extreme observed cases: 
(a) flat X-ray spectrum, broad optical emission-line type 1 
AGN continua, and (b) steep X-ray spectrum, narrow 
optical emission-line type 1 AGN continua. If a peculiar 
class of AGNs does exist, with intrinsic UV to X-ray contin- 
uum switching between these two extreme observed cases 
(as suggested by the behavior observed in 1H 0419 — 577, 
Guainazzi et al. 1999, Turner et al. 1999; NGC 4051, 
Guainazzi et al. 1998; RE J2248 — 511, Puchnarewicz et al. 
1995), then sources belonging to this class would show 
emerging spectra of the types (a) and (b) alternately. 

2.3.1. Flat X-Ray Ionizing Continuum 
Typical O-UV to X-ray continua of broad-line type 1 
AGNs (i.e., a £ ~ 1, ot ox ~ 1.3) transmitted by photoionized 
gas show deep O vn-O vm edges, since those are usually the 
most abundant oxygen ions. In these conditions, when the 
oxygen is mainly distributed in highly ionized states, carbon 
is almost fully stripped and neon is mostly distributed 
between Ne ix and Ne xi. Conversely, high-ionization iron 
ions are not much populated, the main ions being Fe xiv- 
Fe xx (Nicastro et al. 1999). In Figure 3, we show 0. 1-5 keV 
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Energy (keV) 

3 0.1-5 keV spectra emerging from an outflowing wind photoionized by a typical flat X-ray, broad-line type 1 AGN continuum, for log U — 0.5, 1, 

l 5 v _ 1000 km s ~\o v = 500 kms" 1 , and N H = 6 x 10 21 cm“ 2 . The low energy turnover is due to absorption from a column of neutral gas, to mimic 

typical Galactic absorption. 


spectra emerging from an outflowing wind photoionized by 
a typical flat X-ray spectrum, broad optical emission-line 
type 1 AGN continuum (Ferland 1997), for three different 
values of log U = 0.5, 1, 1.5. We used an outflowing velocity 
of v oul = 1000 km s* 1 , which is typical for the UV out- 
flowing component in warm absorbers (Mathur et al. 1995) 
and which corresponds to z = ~ 0.003. For the dispersion 
velocity, we adopted a value o v = 500 km s 1 . The column 
density of photoionized gas along the line of sight is fixed at 
JV h = 6x 10 21 cm -2 , which is the mean value in the Rey- 
nolds (1997) sample. We also added absorption from cold 
gas of column density A^ old = 3 x 10 20 cm" 2 to simulate a 
typical Galactic absorption along the line of sight. In each 
panel of Figure 3 we plot the transmitted spectra corre- 
sponding to a given value of log U . 

Moving from log U = 0.5 to log U = 1.5 (i.e., from low 
to high ionization. Fig. 3, top to bottom) the O vii K edge 
opacity decreases, due to the decreasing of the O vn/O vm 
relative abundances ratio. Conversely the optical depth of 
the O vm K edges increases until the O vii — ► O vm photo- 


ionization rate becomes lower than the O vm -* O ix one. 
In Figure 3 (bottom) (log U = 1.5) only the O vm K edge is 
visible, whereas the gas is almost completely transparent at 
the O vii K edge energy. In the lower ionization case 
(log U = 0.5; Fig. 3, top\ a deep O vii K edge is instead 
present, along with the 0.49 keV K edge of C vi. Ne ix K 
edge is marginally visible because of the low column density 
adopted. 

Strong O vii Ka and K/f and O vm Ka resonance absorp- 
tion lines at 0.574, 0.653, and 0.665 keV, respectively, as well 
as C v, C vi, Ne ix, and Ne x Ka and K/? lines, are present in 
these spectra with intensities and equivalent widths depend- 
ing weakly on the ionization state. On the other hand, 
intensities and equivalent widths strongly depend on the 
dynamics of the absorbing clouds of gas along the line of 
sight and on its column density. For the case under con- 
sideration (a v = 500 km s" 1 , N h = 6 x 10" 1 cm 2 ) the 
equivalent widths of Ka and K/? He-like and H-like carbon, 
oxygen, and neon are in the range 2-8 eV. In the medium 
and higher ionization cases, Ka and KjS lines from Mg and 
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Si are also present, along with two complex systems of 
S xi-S xv and Fe xv-Fe xx L lines around 0.3 and 1 keV, 
respectively. Their equivalent widths in the higher ioniza- 
tion case are of the order of 1 eV for each line. We note that 
the energies of the strong Ka and K/? Ne ix and Ne x lines, 
as well as those of a number of Fe xvii and Fe xvm lines, are 
all close to each other and to the energy of the O vm K 
edge. Caution should then be used when fitting low- and 
medium-resolution spectra of warm absorbers with a 
phenomenological multiedge model. Neglecting the pres- 
ence of this system of absorption lines (which may reach 
total equivalent width of up to a few tens of eV for a v = 
1000 km s" 1 ) and leaving the optical depths of each edge 
free to vary independently in these fits may result in good 
fits from the / 2 point of view, but in overestimating the 
O viii optical depths. 

2.3.2. Steep X-Ray Ionizing Continuum 

When the ionizing continuum is steep below 2 keV (ce £ > 
2) and flattens above this energy to a £ ~ 0.7-1. 5, as in many 
NLSyls, the resulting ionization structure can be quite dif- 
ferent from the ones described in the previous section. In 
fact, in these condition C, O, and Ne would be almost fully 
stripped in the gas, whereas iron would still be distributed 


in medium-low ionization states, giving rise to emerging 
spectra with no sharp absorption features below - 1 keV 
and a broad absorption structure between 1 and 2 keV due 
to a complex mixture of L resonance absorption lines and 
edges from Fe x-Fe xx. 

As an example, Figure 4 shows two spectra emerging 
from an outflowing wind photoionized by a typical NLSyl 
continuum (see Leighly et al. 1997), consisting of a black- 
body with kT = 0.13 keV and a power law with a £ = 0.9. 
The relative normalization between the soft and the hard 
components is such that the ratio between the bolometric 
luminosity of the blackbody component and the 2-10 keV 
luminosity is [L B o 1 (BB)/L 2 _ 10 (PL)] ~ 10 (a mean value in 
the Leighly et al. sample). We assume also in this case 
t out = 1000 km s* 1 , cr v = 500 km s“ l , and a column density 
of N% old = 3 x 10 20 cm~ 2 neutral gas to mimic a typical 
Galactic absorption. In Figure 4 (top) we plot the spectrum 
emerging from gas with log U = -0.5 and log JV H = 22. 
Both absorption lines and edges form the H-like and 
He-like ions of C, O, and Ne, along with a number of Si L 
lines around 0.3 keV, are imprinted on spectra emerging 
from photoionized gas with this (or lower) value of U , 
similar to what is found in the case of flat X-ray ionizing 
continuum; see Figure 3. Figure 4 ( bottom ) illustrates the 
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Fig. 4.— Spectra emerging from an outflowing wind photoionized by a typical NLSyl continuum and with v out = 1000 km s’ 1 and a v = 500 km s ‘ l . The 
ionizing continuum used consists of a blackbody with kT = 0.13 keV and a power law with a E = 0.9. The ionization parameter and gas column density are 
log V = —0.5, log iV H = 22 and log U = 1, log N H — 23 in the upper and lower panels, respectively. 
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case of a much higher ionization (log U = 1). Here we use a 
column density an order of magnitude higher than in the 
previous case to highlight the peculiarities of this model : no 
deep absorption edge is visible in this spectrum, despite the 
very high column density of N H = 10 23 cm' 2 , whereas a 
complex system of strong absorption lines between 1 and 2 
keV is clearly visible, due to almost 60 L lines of midionized 
iron and Ka and K/? lines of Mg, Si, and S. The iron L lines 
around 1 keV are all very close in energy, and their equiva- 
lent widths span from a few eV for a v = 500 km s 1 up to 
tens of eV for a v = 1000 km s _1 . The total equivalent width 
due to these lines between 1 and 1.5 keV can be as high as a 
few hundred eV. 

Let us now suppose that the ionized gas in NLSyl has 
physical properties similar to the warm absorbers in flat 
X-ray spectrum type 1 AGNs, i.e., similar electron density 
and similar distance to the source of the ionizing radiation 
in units of gravitational radii. Pounds et al. (1995) and Laor 
et al. (1997) suggested that steep X-ray spectrum, narrow- 
line quasars are emitting close to the Eddington luminosity, 
at a rate, say, a factor of 10 higher than flat X-ray spectrum 
broad-line quasars. This means that the black hole mass of 
steep X-ray quasars is a factor of 10 smaller than that of flat 
X-ray quasars with a similar luminosity (in agreement with 
variability studies, Fiore et al. 1998a, and with detailed 
modeling of the optical UV X-ray SED of NLSyls, Puch- 
narewicz, Mason, & Siemiginowska 1998). The remarkable 
uniformity of ionization parameters found by Reynolds 
(1997) and George et al. (1998) in their samples of Seyfert 1 
galaxies and quasars of luminosity spanning 3 orders of 
magnitude implies that the product R 2 n e scales linearly with 
the luminosity in these objects. If this correlation holds for 
NLSyls too, then we can estimate the R 2 n e factor of a 
NLSyl of a given luminosity using, for example, the Rey- 
nolds estimate but scaling the distance between the 
absorber and the ionizing source by the above factor of 10. 
For a NLSyl of luminosity 10 43 ergs s _1 we find R 2 n e ~ 
10 40 cm' \ using the typical ionization parameter found by 
Reynolds (1997). Using this value and the correct ionizing 
continua we expect for a typical NLSyl a log U of 0.8- 1.2 
and therefore a spectrum similar to that plotted in Figure 4 
(bottom). 

3. DATA AND SIMULATIONS 

To compare our models with the available data and to 
make predictions on what we will be able to observe with 
the next generation of X-ray spectrometers (i.e., gratings 
and calorimeters), we folded our photoelectric plus resonant 
absorption model with the responses of the ASCA SIS, 
AX AY ACIS/HETG, and a high-resolution, high effective 
area calorimeter as the baseline for the future mission 
Constellation-X. The simulations were performed using 
model parameters as deduced from the ASCA data of the 
Seyfert 1 galaxy NGC 985 and the NLSyl IRAS 
13224—3809. We discuss these two cases in turn. 

3.1. Flat Ionizing Continuum : The Case of NGC 985 

NGC 985 is a Seyfert 1 galaxy known to host an ionized 
absorber (Brandt et al. 1994). Nicastro et al. (1998) observed 
the source with ASCA and confirmed the presence of an 
ionized absorber along the line of sight. They find that a 
single-zone photoionization model (accounting for photo- 
electric absorption only, model a) is not able to explain 
completely the complex absorption structure seen in the 


ASCA spectrum (the continuum has been parameterized as 
a power law with a £ = 1.3 up to 2.3 keV, whereas above this 
energy the spectrum flattens by Aa E = 0.6, the ionized 
absorber parameters are log U = 0.93 ± 0.09, log N H = 
21.95 ± 0.05, and z abs = 0.0232 which implies a 
maximum outflowing velocity % 15,000 kms -1 ). Analysis of 
the ASCA SIS residuals shows a significant deficit of counts 
around - 1 keV (Fig. 5a). The equivalent width of this 
feature (estimated by adding a Gaussian absorption line to 
the continuum model) is — (1 3 ± 4.0) eV. 

We computed the spectrum emerging from gas photoion- 
ized by the NGC 985 optical to X-ray continuum using our 
photoelectric plus resonant absorption model. We extrapo- 
lated the soft X-ray power law, as measured from the ASCA 
data, through optical wavelengths. This is consistent with 
both the measured O-UV spectral index a 0 _ uv =1.4 (de 
Vaucouleurs et al. 1991) and a ox = 1.4. We adopted for log 
U and log N H the best-fit values given above and an out- 
flowing velocity of 1000 km s”\ which is the mean mea- 
sured value for the known O-UV/X-ray ionized absorbers 
(Mathur et al. 1994, 1995, 1997, 1998) and is consistent with 
the estimate of the warm absorber redshift in this source. 
We used a dispersion velocity of 1000 km s' l . The model is 
plotted in Figure 5c. It shows strong resonance absorption 
lines between 0.85 and 1.2 keV due to the redshifted Ka and 
KjS Ne ix and Ka Ne x lines. The blend of these three lines 
has a total equivalent with of about 20 eV, similar to the 
equivalent width of the line feature in Figure 5a. We folded 
this model with the v4SC4-SIS response and simulated a 40 
ks observation of NGC 985. We fitted the resulting spec- 
trum with model a. The ratio between the simulation and 
best-fit model a is shown in Figure 5 b; it resembles strik- 
ingly the residuals to the “real” spectrum in Figure 5a. 
Although this is certainly not a self-consistent, quantitative, 
and unique description of the spectrum of NGC 985, it 
indicates that resonant absorption may well be the cause of 
the feature observed at about 1 keV. 

The spectrum in Figure 5c also shows strong absorption 
Ka and K/J lines from O vn and O vm at 0.57 and 0.65 keV. 
There is no evidence for an absorption feature (in addition 
to photoelectric absorption) at these energies in the ASCA 
spectrum of NGC 985 (Fig. 5a). Contribution from emission 
lines due to recombination could be important at these 
energies, and so they could help in reducing the depth of the 
absorption lines. Line emission strongly depends on the 
covering factor, and therefore an accurate estimate of line 
emission intensity can provide information on the geometry 
and distribution of the ionized gas. However, the above 
energy is close to the lower threshold of the SIS instrument, 
where the calibration is most uncertain and problems due 
to the residual dark distribution (RDD) are large, and 
where the CCD energy resolution is 5-6 only. For a quanti- 
tative study of the O vii and O vm emission/absorption 
lines we should await instruments with better energy 
resolution and coverage. 

As a next step we folded the photoelectric plus resonant 
absorption model of Figure 5c with the responses of AXAF 
ACIS/HETG and the baseline Constellation-X calorimeter, 
which has a factor of 2 worse energy resolution but a factor 
of 100 higher effective area (at 1 keV). The simulations are 
shown in Figures 6a, 6b, 6c and 7a, 7b, 7c for three different 
values of the dispersion velocity: c v = 100, 500, and 1000 
km s' 1 for panels (a), (b), and (c), respectively. The first- 
order AXAF ACIS/MEG (medium energy grating) Simula- 
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Fig. 5. — (a) Ratio between the co-added ASCA SISO and SIS1 data of NGC 985 and the best-fitting model a. (b) Ratio between a 40 ks ASCA SIS 
simulation, obtained using the model plotted (bottom; see text for details), and model a. (c) The spectrum emerging from gas photoionized by the NGC 985 
SED, when accounting for both photoelectric and resonant absorption. We adopted outflowing and dispersion velocities of 1000 km s “ l . 


tions (Fig. 6) were performed using the MARX simulator 
(Wise et al. 1997, MARX User Guide, version 1.0 8 ). The 
exposure time is 80 ks. To simulate the calorimeter spec- 
trum (Fig. 7), we wrote a simple program to fold our models 
with the response matrix 9 and add statistical noise. The 
exposure time in the Constellation-X calorimeter simula- 
tions is 10 ks. 

Resonance absorption lines of width >200 km s -1 are 
detected in the AXAF ACIS/HETG first-order spectrum. 
All the strongest lines in Figure 5c are resolved in the 
spectra of Figures 6b and 6c. AXAF HETG spectra of 
Seyfert 1 galaxies will then establish whether resonant 
absorption is important in the ionizing gas and then 
provide estimates of the gas dispersion velocity and bulk 
motion. This will allow unambiguous solution of the open 

8 http ://asc.harvard.edu/prop_port.html # marker_marx. 

9 ftp://legacy.gsfc.nasa.gov/htxs/calor_w _grat.rsp.gz. 


and controversial issue on the identification between the 
X-ray and the UV ionized absorbers. However, the signal to 
noise will be poor in most cases: the O vii Ka/K/? and Ne ix 
Ka/KjS line ratios, which can be used for plasma diagnostic, 
cannot be estimated with uncertainty <100%; line widths 
and profiles that give information on the geometrical and 
dynamical configuration of the absorber can be roughly 
measured but not studied in detail. The HETG can obtain 
spectra of comparable quality for about 100 AGNs. 

The large effective area of the Constellation-X calorime- 
ter along with a good energy resolution will permit attain- 
ment of high signal-to-noise spectra, of a quality similar to 
optical spectra of AGNs, of sources of the brightness of 
NGC 985. This will allow application of plasma diagnostics, 
accurate measurement of both the chemical composition 
and ionic distribution of the gas, discrimination between 
collisional ionization and photoionization, and study of the 
time-evolving behavior of the relevant physical quantities 
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Fig. 6 — Simulations of 80 ks 0.4- 1.2 keV first-order AXAF ACIS/MEG spectra of NGC 985. A dispersion velocity of<r„ = 100, 500, and 1000 km s' 1 has 
been used in (a), ( b ), and (c), respectively. 


(temperature, relative ionic abundances, etc.), which in turn 
will allow constraint of important physical and geometrical 
parameters (electron density and radial distance; Nicastro 
et al. 1999). The Constellation-X calorimeter could provide 
moderate signal-to-noise spectra, similar to those in Figure 
6, of sources down to a flux limit of a few 10“ 13 ergs cm 2 
s _1 , which increases the sample of possible targets to a few 
tens of thousand sources. AGNs with luminosity in the 
range 10 43 -10 44 ergs s" 1 can be studied up to z = 0. 3-0.5, 
whereas more luminous (10 46 -10 47 ergs s “ l ) quasars will be 
available up to z - 3 (the limit is given by the bandwidth of 
the Constellation-X calorimeter more than by its sensitivity, 
since features at 0.5—1 keV are redshifted to 0.12-0.25 keV 
for z = 3). This will permit the study of the evolution of 
physical state, chemical composition, and geometry of the 
absorbing gas with redshift and luminosity. Statistical 
studies on the distribution of the absorber density, ioniza- 
tion state, chemical composition, and distance from the ion- 
izing source will be possible, which can give information on 
the absorber origin. 


3.2. Steep Ionizing Continuum: The Case of IRAS 
13224-3809 

We chose to test our model on the NLSyl galaxy IRAS 
13224 — 3809, which showed in an .4SC/1 observation a 
deep absorption feature in the interval 1-2 keV, interpreted 
by Leighly et al. (1997) in terms of resonant absorption from 
gas outflowing with v/c ~ 0.5 

Figure 8a shows the ratio between the co-added ASCA 
SISO and SIS1 data and the best-fitting model composed 
by a blackbody plus a power law, with best-fit param- 
eters given by kT = 0.134 keV, <x E = 1.16, and 
[L Bo 1 (BB)/L 2 _ 10 (PL)] = 9.8 (fully consistent with both the 
Leighly et al. 1997 and the Hayashida 1997 estimates). A 
deep and broad deficit of counts is evident between 0.9 and 
2 keV (also see Fig. 1 in Leighly et al. 1997). The absorption 
feature equivalent width (estimated by adding a Gaussian 
absorption line to the continuum model) is large: 
EW « -(180 ± 30) eV. We note that, even if the deficit of 
counts occurs approximately where the blackbody and the 
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Fig. 7.— Simulations of 10 ks spectra of NGC 985 taken with the Constellation-X calorimeter. A dispersion velocity of a v = 100, 500, and 1000 km s' 1 
has been used in (a), (f>), and (c), respectively. The cutoff at low energy is partly due to the inclusion at low energy of neutral absorption from the Galactic 
column along the line of sight (JV H = 3 x 10 20 cm ” 2 ). 


power-law components match each other, we were not able 
to model the absorption-like feature without invoking a 
third component. 

We computed the spectrum emerging from gas photoion- 
ized by the IRAS 13224—3809 UV to X-ray continuum for 
a number of values of log U (from -0.5 to 1.5), log N„ (from 
10 22 to 10 24 cm -2 ), and a v (from 100 to 1000 km s - '). We 
extrapolated down to 13.6 eV and up to 100 keV the best-fit 
continuum to the ASCA data of IRAS 13224—3809. This is 
in good agreement with the optical UV to X-ray measures 
of this source, which suggests that the thermal disk com- 
ponent in this object could be present at higher energies 
than usually found in soft X-ray flat, broad-line type 1 
AGNs (Mas-Hesse et al. 1994). We adopted for the hard 
X-ray power law the best-fit spectral index <x E = 1.16. This 
value is in excellent agreement with the mean estimate of 
Hayashida (1997). However, as Brandt et al. (1997) pointed 
out, the spectral index of the hard power law of IRAS 
13224—3809 is highly variable, ranging between 0.7 and 
1.3. The value that we used in our model represents there- 


fore only a mean value. We stress however that the exact 
shape of the hard power law does not affect strongly the 
distribution of the relative ionic abundances in the gas for 
elements fighter than neon, which is instead mainly driven 
by the presence of a strong soft excess below ~2 keV. Our 
results are therefore rather insensitive to the exact slope of 
the hard power law. As in the previous case we adopted an 
outflowing velocity of 1000 kms' 1 . We folded these models 
with the ASCA-SIS response and then fitted them with a 
blackbody plus power-law model. As an example Figure 8 b 
shows the ratio between an 80 ks simulation obtained using 
log N h = 23.5, log U = 1, outflowing and dispersion 
velocities of 1000 km s -1 (model shown in Fig. 8c), and the 
best-fitting blackbody plus power-law model. We note that 
all the relevant absorption fines lie in a spectral region in 
which the intrinsic continuum is still dominated by the soft 
blackbody component (Fig. 8c, but see also Fig. 4). The 
residuals resemble those in Figure 8 a. The large deficit of 
counts between 0.9 and 2 keV is due to the presence of a 
large number of iron L lines. The co-added equivalenl width 
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Fig. 8. — (a) Ratio between the co-added ASCA SISO and SIS1 data of IRAS 13224-3809 and the best-fitting blackbody plus power-law model, (b) Ratio 
between a 80 ks ASCA SIS simulation, obtained using the model plotted ( bottom ; see text for details) and the best-fitting blackbody plus power-law model, 
(c) The spectrum emerging from gas photoionized by the IRAS 1 3224 — 3809 SED, when accounting for both photoelectric and resonant absorption. We 
adopted outflowing and dispersion velocities of 1000 km $“*. 


of these lines is about 130 eV, similar to the equivalent 
width of the absorption feature in Figure 8a. 

As for the case of NGC 985, we plot in Figures 9a, 9 b, 9c 
and 10a, 10b, 10c three simulations obtained convolving the 
spectrum in Figure 8c with the AX4F/MEG and 
Constellation-X calorimeter responses. The simulations are 
shown for three values of the dispersion velocity, a v =100, 
500, and 1000 km s -1 (Figs. 10a, 10b, and 10c, respectively). 
Integration time was 100 ks for the first-order AX AF / MEG 
spectrum and 20 ks for the Constellation-X calorimeter 
spectrum. 

Many single Fe L lines are well resolved even in spectra 
with dispersion velocity of 1000 km s -1 . AXAFj HETG 
spectra of low-redshift NLSyls with luminosity L = 10 43 - 
10 44 ergs s -t will permit unambiguous establishment of 
whether an ionized absorber is indeed responsible for the 1 
keV deficit of counts seen in many moderate-resolution 
X-ray spectra of these sources. The low signal-to-noise ratio 
of AXAF/HETG spectra, however, will not allow applica- 


tion of plasma diagnostics to infer the physical properties of 
the absorbers. Nevertheless the outflowing and dispersion 
velocities of the gas will be accurately derived, as well as the 
relative Ne to O and Fe to Si and Mg abundances, which 
will put important constraints on the geometry and the 
chemical composition in the AGN environment. The high 
effective area of the Constellation-X calorimeter will allow a 
great increase of the sample of narrow-line type 1 AGNs 
observed at high resolution. The complex iron L edges and 
line absorption system could be detected in high-luminosity 
narrow-line quasars up to redshift <2, or even higher, 
depending on the actual response below the carbon edge 
(0.283 keV), whereas Mg and Si Ka and K/J resonance 
absorption lines could be detected up to z ~ 3. The good 
energy resolution (3 eV), will permit separation of the 
strongest (oscillator strength greater than 0.05) single reso- 
nance absorption lines of the 1 keV iron L edge plus line 
complex, thus allowing for a clear diagnostic of the mecha- 
nism responsible for the ionization of the gas. In particular. 
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Fig. 9.— Simulations of 100 ks 0.4-2 keV first-order AXAF ACIS/MEG spectra of IRAS 13224-3809. A dispersion velocity of er p = 100, 500, and 1000 
km s 1 has been used in (a), (6), and (c), respectively. 


an accurate measure of the relative abundances of the very 
stable Fe xvn Li-like ion could allow clear distinction 
between gas in collisional or photoionization equilibrium 
(Nicastro et al. 1999). Time-evolving ionization can also be 
investigated, thus giving constraints on the electron density 
and the radial distance of the absorbers (Nicastro et al. 
1999). 

4. CONCLUSIONS 

We have calculated spectra emerging from ionized gas 
taking into account both photoelectric and resonant 
absorption. We find that the emerging spectra are strongly 
modified by resonant absorption if the dispersion velocity 
in the gas is >100 km s~L This, or even higher values, is 
not unrealistic, since many broad-line AGNs show absorp- 
tion lines in the UV blueshifted by 200-2000 kms" 1 and 
broadened by similar dispersion velocities. 

In particular we find the following. 

1. The distribution of relative ionic abundances in pho- 
toionized gas dramatically depends on the exact shape of 
the ionizing continuum. 


2. Strong Ka and Kj? resonance absorption lines from 
He-like and H-like ions of C, O, and Ne, along with deep 
C vi, O vn-O viii, Ne ix-Ne x K edges, and Fe xv-Fe xvii 
L edges, are expected to be present in the 0.1-3 keV spectra 
of flat X-ray spectrum quasars emerging from photoionized 
gas along the line of sight. When fitting low- and medium- 
resolution X-ray data with phenomenological multiedge 
models, these resonance absorption lines may produce the 
effect of overestimating the optical depth of the O viii K 
edge. 

3. Conversely, no strong absorption edge is expected to 
be imprinted on spectra emerging from gas illuminated by 
steep ionizing continua, with physical and geometrical 
properties similar to those found in broad-line type 1 
AGNs. However, a large number of strong Fe L, Mg, Si, 
and S Ka and K/? resonance absorption lines (almost 70) are 
predicted between 1 and 2 keV. Fitting low- and medium- 
energy resolution data with continuum models may then 
result in deep smooth negative residuals at these energies. 
We then suggest that the ^ 1 keV absorption feature 
observed by ASCA in several narrow-line Seyfert 1 galaxies 
can be explained in terms of resonance absorption lines. 
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Fig. 10 .— Simulations of 20 ks 0.3-2 keV spectra of IRAS 13224-2809 taken with the Constellation-X calorimeter. A dispersion velocity of <r„ = 100. 500, 
and 1000 km s * 1 has been used in (a), (b), and (c), respectively. The cutoff at low energy is partly due to the inclusion of neutral absorption from the Galactic 
column along the line of sight (A' H = 4.8 x 10 2 ° cm 2 ). 


without requiring relativistic outflowing velocities of the 
gas. 

4. The equivalent width of the predicted resonant 
absorption lines strongly depends on the column density of 
the absorbing gas and its dispersion velocity. Assuming the 
range of N H typically observed and dispersion velocities of 


0-1000 km s' 1 , we estimate equivalent widths spanning 
from a few eV to a few hundred eV. 

F. F. acknowledges support from NASA grants NAG5- 
2476 and NAG5-3039. 
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